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Unsteady  Forces  Acting  on  a  Flat  Plate  Airfoil 
Undergoing  Chord  Deformations 

John  H.  Russell*  and  Chuen-Yen  Chow** 

Abstract 

This  paper  investigates  the  forces  generated  by  a  flat  plate  airfoil  that  undergoes 
prescribed  chord  changes  in  an  otherwise  uniform  flow.  An  inviscid  and  incompressible 
fluid  is  assumed  and  conformal  mapping  techniques  are  used  to  transform  a  deforming 
circular  cylinder  in  the  computational  plane  into  a  flat  plate  airfoil  undergoing  chord 
deformations  in  the  physical  plane.  The  resultant  unsteady  lift  is  determined  and 
discussed. 

I.  Introduction 

The  utilization  of  unsteady  flow  fields  for  lift  enhancement  continues  to  receive  the 
attention  of  both  experimental  and  theoretical  investigators  (Ref  1).  A  significant 
number  of  these  studies  deal  with  the  generation  and  control  of  large  vortex  structures 
in  the  vicinity  of  the  wing  upper  surface.  The  large  lift  increments  achieved  as  a  result  of 
these  vortex  structures  have  been  documented  by  both  experimental  (Ref.  2)  and 
theoretical  (Ref.  3)  efforts.  The  current  study  differs  from  the  aforementioned  in  that 
unsteady  effects  enter  the  problem  only  through  time  dependent  shape  (chord)  change 
undergone  by  a  flat  plate  airfoil  and  through  the  effects  of  the  shed  vortex  wake.  An 
earlier  investigation  (Ref  4)  determined  that  a  degree  of  unsteady  lift  enhancement 
arose  from  prescribed  deformations  by  an  elliptical  cylinder.  That  analysis  showed  this 
increment  to  occur  in  the  absence  of  unsteady  circulation  about  the  body  and  without 
the  associated  shed  wake.  The  current  study  considers  the  circulatory  features  of  the 
flow  that  result  from  continuous  enforcement  of  the  Kutta  condition  at  the  trailing  edge 
of  a  deforming  flat  plate  airfoil,  and  the  effects  that  these  unsteady  circulatory  features 
have  on  the  generated  lift. 


•Major,  USAF,  Instructor,  Dept,  of  Aeronautics,  USAF  Academy,  CO. 

••Professor,  Dept,  of  Aerospace  Engineering  Sciences,  Univ.  of  Colo. ,  Boulder,  CO. 
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II,  The  Computational  Model 

In  their  earlier  work,  Russell  and  Chow  (Ref.  4)  developed  a  computational  model  for 
a  circular  cylinder  deforming  radially  in  a  uniform  flow  that  is  inclined  at  an  angle  of 
attack  a.  For  the  model  shown  in  Figure  1 ,  the  complex  potential  is  given  by 


w(i)  =  U  |<  —  e‘“ 


da 

+  a  —  In  C 
dt 


(1) 


Figure  1.  A  radially  deforming  circular  cylinder  in  a 
uniform  flow  at  an  angle  of  attack. 

Now  the  model  is  extended  to  analyze  the  unsteady  forces  arising  from  prescribed 
chord  deformations  of  a  flat  plate  airfoil.  As  a  result  of  continuous  enforcement  of  the 
Kutta  condition,  the  circulation  about  the  airfoil  is  itself  time  dependent.  The  constancy 
of  circulation,  in  a  global  sense,  thus  requires  a  shed  vortex  wake  in  response  to  the 
unsteady  bound  vortex.  The  added  implications  of  the  unsteady  circulation  are  discussed 
here  and  the  circulatory  terms  are  added  to  the  previously  developed  model. 
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The  corresponding  flow  condition  is  shown  in  Fig.  2b  where,  in  fact,  an  infinite  velocity 
exists  at  the  sharp  trailing  edge.  In  order  to  re-establish  tangent  flow  at  the  trailing 
edge,  the  bound  circulation  must  increase  to  a  value  of  =  To  +  ff .  In  order  for  the  total 
circulation  in  the  flow  to  remain  constant,  vorticity  of  strength,  -K,  must  be  shed  as 
shown  in  Fig.  2c. 

In  as  much  as  the  deformations  are  continuous  in  time,  so  would  be  the  changes  in 
bound  vorticity.  Likewise,  the  shed  vorticity  would  be  continuously  distributed 
throughout  the  wake.  An  analytical  solution  for  this  case  would  be  extremely  difficult, 
hence  a  simplified  model  is  desired.  If  the  time  increment  of  deformation  At  is 
sufficiently  small,  enforcement  of  the  Kutta  condition  can  be  made  arbitrarily  close  to 
continuous.  If,  after  each  small  time  increment,  the  condition  of  tangency  of  flow 
leaving  the  trailing  edge  is  applied  to  arrive  at  a  new  bound  circulation  value,  then  a 
discrete  point  vortex  will  be  shed  into  the  wake  in  a  process  similar  to  that  modeled  by 
Chow  and  Huang  (Ref.  3).  As  the  deformations  continue  in  time,  the  above  sequence  is 
continuously  repeated  thus  leaving  a  series  of  point  vortices  trailing  the  airfoil.  Figure  3 
depicts  the  previously  defined  computational  model,  now  with  the  addition  of  the 
circulatory  features  just  described. 

In  this  analysis  the  flat  plate  is  considered  as  a  limiting  case  of  a  symmetrical  airfoil 
whose  thickness  approaches  zero.  This  avoids  the  difficulty  attendant  with  infinite 
velocities  that  arise  at  the  leading  edge  of  the  flat  plate.  Under  the  Joukowski 
transformation 


2  =  <  + 


< 


(2) 


where 


c(t)  =  a  (t)  —  E 
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i,=  i,  -t 


a  circle  of  radius  a  in  the  C'  plane  transforms  into  a  symmetrical  airfoil  in  the  z  plane 


Figure  3.  Joukowski  transformation  of  a  deforming  circular  cylinder 
into  a  nearly  flat  plate  airfoil  at  an  angle  of  attack. 

If  e  isthought  of  as  a  thickness  parameter  then  the  airfoil  thickness  will  be  0(c)  ,the 
leading  edge  radius  is  0(t2) ,  and  the  chord  is  4c  +  0(t2)  (Ref.  6).  The  point  <  =  c  in  the 
computational  plane  maps  into  z=2c  in  the  physical  plane  which  corresponds  to  the 
airfoil  trailing  edge.  Figure  4  provides  a  more  detailed  view  of  the  area  surrounding  the 
origins  of  the  <'  and  <  coordinate  systems.  The  limiting  process  c  -*0  recovers  the  flat 
plate  airfoil  with  the  leading  edge  of  zero  thickness  located  at  z  =  -2c  in  the  physical 
plane.  Deformations  to  the  airfoil  chord  arise  entirely  through  radial  deformations  of 
the  circular  cylinder.  Under  the  transformation  of  Eq.  (2),  a  freestream  of  velocity  Ua>  in 
the  computational  plane  maps  into  the  same  flow  in  the  physical  plane.  The  circulatory 
features  of  the  inviscid,  incompressible  flow  are  additive  to  those  elements  developed  in 
Ref.  4.  With  the  aid  of  the  Circle  Theorem  (Ref.  7),  the  resulting  complex  potential  is 
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written  in  the  fixed  C  plane 


=  —  c‘“ 


ir  •  da 

•h  —  In  C  +  a  —  In  C  -f- 

2n  dt 


iff  f  ] 

s 

where  the  overbar  denotes  the  complex  conjugate.  The  first  and  third  terms  . 
respectively,  those  for  the  rigid  circular  cylinder  in  a  uniform  flow  and  the  deform 
cylinder  as  developed  in  Ref.  4.  The  second  and  fourth  terms  represent  the  complex 
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potential  for  the  circulation  bound  to  the  airfoil,  ,  as  well  as  the  potential  for  the 
discrete  shed  vortices,  K  ,  and  their  images  within  the  circle.  As  was  previously 
accomplished  for  the  elliptical  cylinder  (Ref.  4),  the  forces  produced  by  the  deformations 
and  the  surrounding  flow  field  are  determined  by  the  unsteady  Blasius  Theorem  as 
developed  in  Ref.  8.  The  appropriate  form  for  this  analysis  is  given  by 


ip  [  /dw\^  d  f  ^ 

X+iV  =  -ip2^—  +  -  {—  +  T  P-JT 


dt 


c  .  \  dw  \  I A 

-  +  ip  Z - 


dt  I  dt 


X  and  Vare  the  resultant  forces  aligned  with  the  x  andy  axes  of  the  physical  z  plane.  F^ 
is  the  vorticity  bound  to  the  airfoil  and  the  term  Sz^  represents  the  first  moment  of  area 
of  the  airfoil  cross-section.  The  lift  and  drag  are  resolved  through  the  transformation 


D  iL  =  (X  +  iY)  e 


In  terms  of  the  computational  coordinate,  Eq.  (A)  is  rewritten  as 


t—dC  -h 
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and 


dz 


(8) 


In  as  much  as  this  is  a  study  of  the  forces  arising  from  prescribed  body 
deformations,  it  should  be  understood  that  the  values  for  a  and  da/dt  as  well  as  C/®  and  a 
are  input  parameters  of  the  problem.  Then  as  posed  by  Eqs.  (6-8),  the  forces  are 
determined  by  analysis  of  the  circulation  elements  of  the  flow  field. 

III.  The  Kutta  Condition 

In  a  reference  frame  fixed  to  the  trailing  edge  of  the  airfoil,  the  complex  velocity 
of  the  flow,  relative  to  the  moving  trailing  edge,  is  given  by 


du>  dw 


dz 

'n 


(9) 


rel  dl  C 

Since,  in  the  computational  plane,  the  trailing  edge  is  moving  to  the  right  at  velocity 


daldt ,  then 


dw  ^  dw  ^  da 

di,  ^rel  dl,  <  =a 


(10) 


where  dw/dC  is  given  by  Eq.  (7).  When  Eqs.  (7),  (8)  and  (10)  are  substituted  into  Eq.  (9) 
and  the  result  is  evaluated  at  +he  trailing  edge,  we  get 


■  ^  miK 

—  2iU  sin  a  + - +  ^  — - 

*  2na  2n 


/ 


a-C 


<-« 


(11) 


The  resulting  infinity  at<  =  cisof  course  the  constraining  factor  on  the  bound  circulation 
strength  and  the  individual  shed  vortex  strengths.  An  indeterminate  form  exists  if  the 
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numerator  of  Eq.  (1 1)  can  be  set  identically  to  zero  such  that 


ir  ^ 

-2iU  sina+ - +  7  - 

"  2na  2n 


2 

a 

a  -  — 


For  a  closed  contour  that  encloses  the  airfoil  and  the  entire  wake,  Kelvin's 
Circulation  Theorem  requires  that,  within  the  contour,  the  circulation  remain  constant. 
Taking  this  constant  value  to  be  that  of  the  circulation  required  of  an  airfoil  of  constant 
chord  4a|j  operating  under  the  same  flight  conditions,  then 


while  the  bound  circulation  T.  is  found  from  the  constraint 


r.=r  -h  y  K.  (1^) 

0  a  ^  j 
j=l 

But  the  circulation  strength  T  of  Eq.  (12)  includes  both  the  bound  vortex,  r^,  and  that 
portion  of  the  shed  vortex  image  system  located  at  the  origin  of  the  C  plane.  Thus 


because 


r  =  r  +  y  K. 

a  ^  j 
j=t 

and  the  constraint  of  Eq.  (14),  we  can  say  F  =  Fq  in  Eq.  (12).  Then  the  following 
relationship  identically  satisfies  the  Kutta  condition  and  fixes  the  strength  K„^  of  the  mtf’ 
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vortex  shed. 


tm  —  l 

4nU ^(a—a^  sina  —  ^  Kj 


J  J 


j=l  Ua-C,^)(a-i..) 


2  7  7 

a  C, 
m  m 


(a-C  )(a-  O 

ttl  m 


Cm  is  the  position,  in  the  computational  plane,  where  the  vortex  is  introduced  into 
the  flow  field.  In  this  analysis,  the  shed  vorticity  is  assumed  to  originate  in  the  extended 
plane  of  the  trailing  edge.  The  distance  downstream  from  the  trailing  edge  for 
introduction  of  the  vortex  is  a  function  of  the  local  velocity  evaluated  at  the  trailing 
edge.  This  requires  only  determination  of  the  horizontal  velocity  component  for  the  case 
of  the  flat  plate  airfoil.  Due  to  the  indeterminacy  of  Eq.  (11),  the  trailing  edge  velocity  is 
evaluated  alternatively  as 


c^w  cFz  Mc’ 

Another  differentiation  of  Eqs.  (7)  and  (8),  with  evaluation  at  the  trailing  edge,  finds  the 
horizontal  velocity  component  to  be 


c  da  f  c 

U^^=-U  cosa  -h —  /-  — 


c  ^  1 

-Real  —YiK  - — 


n  * 


(a  -  ‘ 


Then  (,'m  is  determined  simply  from  performing  the  inverse  transformation 
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on  the  point 

z  =2c+U.„^t  (19) 

m  i  Ct 

In  practice,  Eq.  (17)  is  evaluated  prior  to  the  introduction  of  the  vortex  and  actually 
represents  the  trailing  velocity  at  the  end  of  the  previous  time  increment.  With  the  shed 
vortex  position  now  determined,  Eq.  (15)  is  used  to  fix  its  strength  within  the  Kutta 
condition  requirement  of  tangent  flow  at  the  trailing  edge.  Note  that  the  strength  of 
this  most  recently  shed  vortex  is  affected  by  the  deformed  geometry  of  the  airfoil 
through  its  new  chord  length  4a.  Additionally,  the  time  dependent  historical  nature  of 
the  unsteady  problem  presents  itself  through  the  effects  of  the  previously  shed  vortices. 

IV.  The  Resultant  Force 

The  lift  produced  by  the  deformations  is  determined  by  solution  of  Eq.  (6).  An 
alternative  form  is  written  by  splitting  this  relationship  into  quasi-steady  and  purely 
unsteady  components.  The  quasi-steady  portion  of  the  Blasius  Theorem  is  written 


The  phrase  quasi-steady  here  only  implies  that  no  explicit  time  dependence  appears  in 
the  given  terms  of  the  Blasius  Theorem.  However,  since  the  complex  potential  has  the 
deformation  term  in  it,  a  degree  of  unsteadiness  is  present  even  for  this  apparently 
steady  term.  Performing  the  indicated  integrations,  one  finds  that 


dc 

dt 


m 

ip  U^r  -2  n p  U ~  \  -  ip  ^  K ,  (u-hiv)  + 


0(c)  (21) 


j  =  i 


The  term  (u  +  iv),^  in  Eq.  (21)  is  the  physical  plane  velocity  induced  at  the  vortex  by  all 
of  the  other  flow  elements  present  and  by  the  singularity  images  within  the  circle.  A 
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general  development  of  the  expression  for  this  induced  velocity  is  included  in  the 
Appendix . 

For  the  obviously  unsteady  contributions  from  Eq.  (4),  we're  left  with 


dr 


dSz 


dt 


+  ip 


(22) 


The  first  term  on  the  right,  the  time  rate  of  change  of  bound  circulation,  requires  some 
clarification.  If  the  approximation 


dr  V  r  -r  , 

m  m  —  l 


ai  \ 

x) ' 


(23) 


At 


is  used,  then  for  the  sequence  required  by  the  invariance  of  circulation  in  the  global 
sense 


we  get 


^  r^-r^=-K^ 


ro  =  r2  +  K^  +  K^  ^  r2-r^=-K^-K2 


^2-^1=  -^2 


r  -r  =  ~K 

m  m-l  m 


dr 


ai  \  K 

_ a  \  m 

dt  Jn  At 


(24) 


(25) 


Equation  (25)  is  merely  a  statement  that  the  instantaneous  rate  of  bound  circulation 
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change  is  equal  to  the  strength  of  the  most  recently  shed  vortex  divided  by  the  time 
increment. 

The  second  term  from  Eq.  (22)  is  evaluated  even  for  the  present  case,  since  the  flat 
plate  is  expected  to  be  recovered  only  in  the  limit  e  -►0.  For  non-zero  values  of  e,  the 
first  moment  of  the  area,  Szc.  and  its  time  derivatives  are  expectedly  non-zero  and  are 
evaluated  as  follows. 

If  the  area  is  defined  as  the  surface  integral 


S  = 


dxdy 


then  application  of  Green's  Theorem  in  its  complex  form  (Ref.9) 


j  j  ^  cixdj'  =  -  j  I  F{z,  z)dz  (26) 

results  in,  for  the  profile  cross-sectional  area 

S=-M  ~z  dz  (27) 

2  Je 

C  is  now  the  closed  contour  defining  the  airfoil  perimeter.  Similarly,  for  the  first  moment 
of  the  area  as  defined  by  the  surface  integral 


utilization  of  Eq.  (26)  produces 


Sz 


z  z  dz 


(28) 


,v 


U 


■V*’’ 


13 


^  E 1  n  f 


USAFA-TR-85-11 


For  the  first  moment  of  area,  evaluation  of  Eq.  (28)  gives,  through  first  order  in  e 

Sz  =  -2  nc^t  +0(e^) 
c 

Evaluationof  Eq.  (21)as  e  -*0  results  in 


K 


dc 


m  (*c  ^  a 

X  +iY„  —  ip2c  —  +i8npU  c  —  sma+ip>  K  — 
2  2  ^  M  ^  dt  J  dt 


j  =  i 


c  c 

VI 


+  0(e) 


(29) 


Together,  Eqs.  (21)  and  (29)  constitute  the  solution  for  the  "near"  flat  plate  that  is 
undergoing  chord  deformations.  For  the  limiting  case  where  e-»0,  Eq.  (21)  and  (29)  are 
equivalent  to  the  results  achieved  by  Sarpkaya  (Ref.  10),  who  analyzed  an  inclined  flat 
plate  of  constant  chord  without  bound  circulation.  Sarpkaya  assumed  a  Kutta  condition 
at  both  the  leading  and  trailing  edges  that  resulted  in  finite  velocities  at  both  ends.  In 
effect,  the  present  analysis  assumes  a  finite  velocity  to  occur  at  the  leading  edge  because 
of  its  extremely  small,  but  non-zero,  radius  of  curvature.  Finiteness  of  ailing  edge 
velocity  is  insured  in  the  traditional  manner  through  enforcement  of  the  Kutta  condition. 
Because  the  velocity  remains  finite  at  the  leading  edge,  the  paradox  of  Cisotti  (Ref.  1 1) 
does  not  appear  in  the  current  study. 

For  the  purely  steady  case,  in  which  no  chord  deformations  occur,  no  discrete 
vortices  are  shed  as  per  Eq.  (15).  Thus  only  the  first  steady  term  from  Eq.  (21)  would 
remain  and  the  lift  predicted  by  the  Kutta  -Joukowski  Theorem  is  recovered. 

If  the  lift  and  drag  are  normalized  by  the  steady  state  lift  that  is  generated  by  the 
original  airfoil  geometry,  then  the  relative  effects  of  the  deforming  airfoil  can  be 
determined.  If  this  steady  lift  is  that  derived  from  the  Kutta-Joukowski  Theorem,  then, 
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for  the  original  airfoil  geometry,  it  is  given  by 

where,  from  Eq.  (13) 


Then  the  lift  and  drag  information  derived  from  the  previous  analysis  is  displayed  in  its 
normalized  form 


(30) 

% 

In  this  form,  however,  the  forces  arising  from  the  unsteadiness  of  the  actual 
deformations  as  well  as  the  steady  forces  arising  from  the  newly  deformed  geometry  are 
not  separable.  These  two  contributions  become  decoupled  when  the  computed  lift  and 
drag  are  normalized  by  the  steady  lift  that  would  be  achieved  by  the  airfoil  in  its 
deformed  state.  Now  the  normalization  is  accomplished  in  the  form 

P  +  (31) 

L 

s 

where 

L  =pu  r 

s  OD  s 

and 


r  =  4  aU  a  sin  a 

s  ® 


Effectively,  in  the  case  of  the  flat  plate,  the  only  difference  in  the  two  normalization 
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schemes  is  that  the  second  method  incorporates  the  deformed  chord. 

Computationally,  the  lift  is  determined  at  each  time  step  as  follows.  For  a 
specified  chord  deformation  rate/freestream  velocity  ratio  (EVR)  and  time  step,  the  new 
chord  length  is  determined.  As  discussed  earlier,  due  to  the  newly  deformed  geometry, 
the  original  circulation  must  be  increased  by  an  amount  sufficient  to  maintain  tangential 
flow  at  the  trailing  edge.  A  discrete  vortex  is  shed  at  the  end  of  the  deformation  whose 
position  is  fixed  by  evaluation  of  Eq.  (17)  for  the  trailing  edge  velocity  and  substitution  of 
this  value  into  Eq  (19).  The  most  recently  shed  vortex  strength  is  then  determined  by 
solution  of  Eq.  (15).  The  complex  velocities  of  each  of  the  shed  wake  vortices  are 
computed  as  described  in  the  Appendix  and  used  in  Eqs.  (21)  and  (29)  to  find  the 
instantaneous  resultant  forces.  Through  the  transformation  of  Eq.  (5),  the  lift  resulting 
from  the  deformations  is  determined  and  then  normalized  by  Eqs.  (30)  and  (31).  The 
wake  vortices'  new  positions  in  the  flow  are  determined  as  discussed  in  the  Appendix. 
The  trailing  edge  velocity  is  evaluated,  and  another  prescribed  deformation  occurs.  The 
newest  vortex  is  introduced  into  the  wake  and  the  computational  process  is  repeated. 

V.  Results  and  Conclusions 

Figures  5  and  6  provide  graphical  description  of  the  lift  enhancement  derived  from 
different  chord  deformation  rates.  Figure  5  shows  the  imaginary  component  of  Eq.  (30) 
for  an  angle  of  attack  of  8°  The  curves  shown  here  are  representative  of  those  for  all 
practical  angles  of  attack.  Because  of  the  normalization  process,  the  angle  of  attack 
dependence  of  the  L/Lq  is  effectively  removed  for  angle  of  attack  generally  less  than  15°. 
For  greater  angles  of  attack  at  larger  EVR  values,  the  closer  proximity  of  stronger  discrete 
vortices  in  the  wake  tends  to  give  slight  normalized  lift  variation  with  angle  of  attack. 
The  normalization  of  lift  values  given  in  Fig.  5  shows  that  one  obvious  effect  of  chord 
stretching  is  an  increase  in  generated  lift.  This  manner  of  data  presentation,  however. 
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Time 

aolUmSec 


Figure  5.  Unsteady  L/Lq  ratio  generated  by  a  flat  plate  airfoil 
undergoing  chord  stretching.  a=8° 

combines  those  lift  enhancement  features  resulting  from  both  the  new  geometry  and 
the  unsteadiness  that  produces  the  new  geometry.  The  normalization  of  Eq.  (30)  does 
provide  an  accurate  time  history  of  the  expected  lift  enhancement  that  would  result 
from  chord  stretching.  The  results  cited  in  Ref.  4  indicated  that  larger  EVR  values  resulted 
in  more  favorable  lift  enhancement.  Figure  5  verifies  that  characteristic,  but  in  this  case, 
for  a  body  with  circulation  based  lift  dominating. 

Figure  6  provides  a  view  of  the  purely  unsteady  effects  that  are  derived  from  the 
chord  stretching.  Because  the  normalization  of  Eq.  (31)  also  effectively  removes  the 
angle  of  attack  dependence,  the  curves  shown  for  an  angle  of  attack  of  8°  are 
representative  of  those  for  all  practical  angles  of  attack.  The  earlier  discussion  about 
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greater  angles  of  attack  with  large  EVR  values  still  applies.  The  curves  indicate  that  chord 
stretching  itself  provides  an  unsteady  increment  to  the  lift.  The  more  rapidly  the  chord  is 
made  to  stretch,  the  larger  is  the  increment.  The  time'  behavior  of  this  normalized  lift 
shows  that  the  unsteady  contribution  becomes  a  proportionally  smaller  part  of  the 
overall  lift  generation.  This  results  from  the  fact  that  the  total  lift  developed  by  the 
deformations  is  continuously  normalized  by  a  "steady”  lift  value  that  increases  with 
time. 


Figure  6.  Unsteady  L/Ls  ratio  generated  by  a  flat  plate  airfoil 
undergoing  chord  stretching.  a  =  8° 
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APPENDIX 

THE  MOTION  OF  THE  DISCRETE  VORTEX  WAKE 

This  analysis  has  shown  that  the  forces  arising  from  the  unsteady  deformations 
of  an  airfoil  require  evaluation  of  the  velocities  induced  at  each  of  the  discrete  vortex 
centers.  From  Milne-Thompson  (Ref.  1 2),  the  complex  velocity  induced  at  the 
vortex  by  all  of  the  other  elements  of  the  flow  field  is  given  by 

d 

(u-iv)  =  lim  —  w  ~  i  —  ln(z  —  z  )  (A.1) 

\  dz  [  2n  "  J 

n 

Here,  wfzj  is  the  complex  potential  of  all  of  the  flow  elements  including  the  discrete 
vortices  and  their  associated  image  systems  as  shown  in  Fig.A.1.  Equation  (A.1)  is 
written  alternatively  as 
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Equation  (A.2)  becomes 


^  «l  2..  dt^ 


"*  iK 

I  — 


2 

•  a 
^n-  = 


'n  iK  , 

+  I  —  -r^ 

2n  7' 

■'*'*  ^“S 


'  '  2n  1 7  _  7  rfz  2  —  2 


(A.3) 


Evaluation  of  the  first  and  largest  bracket  term  is  straightforward  for  the  known 
positions  of  the  vortex  centers.  The  second  bracketed  quantity  of  Eq.  (A.3)  is 
analyzed  by  assuming  that  the  term  within  the  brackets  can  be  expanded  in  a  Taylor 
series  about  the  point  <’„  .Thus 


dC,  fd(\  ■  ■  d  \  1  ■■  2 

dz  \  dz  /  ^  U 


which  then  becomes 


dC,'  (d(  \  ■  .  (  d^z  X  /di,' 

dz  ~  (  dz  Y  ‘■'•’I  j/2)  .(  dz 
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(A.4) 
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For  the  term  1l(z-ZrJ,  consider  first  the  expansion 


«-<„)■ 


i—) 


+  0(^ 


(A.5) 


If  the  reciprocal  of  both  sides  of  Eq.  (A.5)  is  taken,  it  becomes 


Applying  the  Binomial  Theorem  to  Eq.  (A.6)  results  in 


+  0((;' 


(A.7) 


If  Eqs.  (A.4)  and  (A.7)  are  combined,  then  the  second  bracketed  term  in  Eq.  (A.4) 
becomes 
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(A.8) 


For  the  Joukowski  transformation 
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the  limiting  value  of  Eq.  (A.8)  is  then  found  to  be 
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(A.9, 


Combining  all  of  the  terms,  the  velocity  induced  at  the  nt^  vortex  becomes 
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The  induced  velocity,  determined  by  solution  of  Eq.  (A.IO),  is  then  used  in  the 
computation  of  the  unsteady  aerodynamic  forces  resulting  from  the  deformations. 

With  time,  the  point  vortices  are  swept  away  from  the  airfoil  as  though  they 
were  fluid  particles.  The  trajectory  of  each  vortex  center  is  determined  by  using  the 
predictor-corrector  integration  scheme  of  Giesing  (Ref.  13).  z(to)  represents  the 
position  of  the  nf^  vortex  in  the  physical  plane  at  time  to,  then  the  predictor  step  of 
the  integration  scheme  places  the  new  position  Zp  at 


=  zit^ 


At 


where  the  overbar  denotes  the  complex  conjugate.  The  velocity  induced  at  this  new 
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predicted  point  is  computed  and  the  final,  corrected  position  at  time  + is  found 
from 


u  [z (r^j)] -  i  u[z (tg)]  +  u  [z  1  -  i  V  [z  ] 


Far  from  the  airfoil,  the  dominant  term  in  the  induced  velocity  relationship 
becomes  the  freestream  flow,  which  remains  constant  with  time,  and  the  effects  of 
the  closest  neighboring  vortices.  Figure  A.2  depicts  the  typical  wake  trajectory  for  the 
case  where  the  prescribed  deformations  result  in  an  increase  in  bound  circulation. 


Figure  A.2.  Discrete  vortex  trajectory  for  typical  defornnation  that 
resulted  in  increase  in  bound  circulation. 

This  case  represents  the  generally  favorable  one  from  a  lift  enhancement  point  of 
view  and  the  wake  is  shown  to  roll  up  as  expected.  It  was  determined  that  for  all 
practical  purposes,  the  wake  itself  provided  little  interaction  with  the  freestream, 
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Unsteady  Surface  Pressure  Measurements  On  A 

Pitching  Airfoil 

John  M.  Walker*  Henry  E.  Helin  **  and  David  C  Chou*** 

Abstract 

Surface  pressure  measurements  were  taken  in  an  experimental  investigation  of 
the  energetic  dynamic  stall  vortices  and  the  associated  unsteady  aerodynamics  generated 
by  and  surrounding  a  6-in.  NACA  0015  airfoil  pitching  at  high  rates  to  large  amplitudes. 
The  airfoil  was  pitched  at  rates  from  230°/sec  to  T380°/sec  and  amplitudes  from  OP  to  60° 
about  its  quarter-chord  axis.  Pitch  rate,  Reynolds  number,  and  the  non-dimensional  pitch 
rate,  a*,  were  varied  to  study  the  effects  on  pressure  and  lift  coefficients.  The 
experimenters  found  that  increases  in  pitch  rate  and  Reynolds  number  had  inverse 
effects  on  the  flow  field  in  the  immediate  vicinity  of  the  airfoil  ;  consequently,  the 
maintenance  of  a  constant  non-dimensional  pitch  rate  produced  very  similar  results. 


I.  Introduction 

As  pointed  out  by  McCroskey,  (Ref.  1)  the  energetic  nature  of  the  unsteady  flowfields 
has  been  a  topic  of  study  for  most  of  the  twentieth  century.  Because  of  their 
complicated,  rapidly  changing  time  dependent  nature  ,  only  recently  has  significant 
progress  been  made  in  theoretical  and  experimental  efforts  to  understand  the  fluid 
mechanics  of  these  flowfields.  A  particular  unsteady  flow  phenomenon,  dynamic  stall 
has  become  increasingly  important  because  of  the  need  not  only  to  predict  its  onset  but 
also  to  understand  its  active  nature  and  exploit  the  associated  increases  in  aerodynamic 
coefficients.  The  importance  of  the  former  has  been  well  known  for  years  with  regard  to 
helicopter  rotors,  turbo-machinery,  and  vertical  axis  wind  turbines.  A  large  part  of  the 
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latter  has  received  impetus  from  Herbst  (Ref.  2)  who  studied  the 
“supermaneuverability"  of  fighter  aircraft.  Understanding  of  highly  unsteady 
flowfields  has  been  hampered  by  two  factors:  1)  the  absence  of  accurate  analytical 
methods  -  short  of  full  Navier-Stokes  solutions  which  are  only  now  becoming  viable  -  to 
predict  these  flowfields  and  2)  a  lack  of  experimental  data  needed  to  set  modeling 
parameters  and  to  verify  theoretical  calculations. 

Most  of  the  experimental  data  on  unsteady  flows  about  pitching  airfoils  to  date 
have  been  restricted  to  sinusoidally  oscillating  models  undergoing  relatively  small 
amplitude  motions  (up  to  ±  fOo)  about  relatively  low  mean  angles  of  attack  (0°  -  15°). 
These  are  typified  by  the  experiments  of  McCroskey  and  Philippe  (Ref.3),  McAlister  and 
Carr  (Ref  4),  Martin,  et.  al,(Ref.  5),  and  Robinson  and  Luttges(Ref.6).  These  data, 
including  flow  visualization,  hot-wire,  and  pressure  measurements,  are  applicable  to 
many  fluid  devices  of  importance.  With  the  exception  of  Robinson  and  Luttges  (Ref.7), 
however,  these  studies  have  not  focused  on  the  fact  that  the  extremely  energetic  nature 
of  the  unsteady  flowfields  could  be  used  to  enhance  performance.  Carr,  et.  al.  (Ref.8) 
have  shown  that  lift,  drag,  and  moment  coefficients  associated  with  unsteady  flowfields 
greatly  exceed  their  static  counterparts.  However,  a  limited  amount  of  experimental 
data  have  been  obtained  from  the  very  large  amplitude  motions  that  will  be  required  if 
the  maneuvers  proposed  by  Herbst  (Ref.2)  using  angles  of  attack  in  excess  of  45°  are  to 
be  implemented. 

Studies  of  moderate  to  targe  amplitude  constant  rate  motions  have  been 
performed  by  Harper  and  Flanigan  (Ref.  9),  who  obtained  force  balance  data  on  a  small 
model  pitching  up  to  30°,  Francis,  et.  al  (Ref  11),  who  obtained  surface  pressure 
measurements  on  airfoils  pitching  up  to  60°,  Deekins  and  Kuebler  (Ref.  1 2),  who  obtained 
flow  visualization  of  dynamic  leading  edge  separation  at  low  Reynolds  numbers  in  a 
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smoke  tunnel,  Daley  (Ref.  13),  who  obtained  leading  edge  dynamic  stall  data  at  low 
Reynolds  numbers,  and  Walker,  et.  al.  (Ref.  14)  and  Helin  and  Walker(Ref.  15),  who 
obtained  flow  visualization  using  a  smoke-wire  and  near  surface  hot-wire  velocity 
magnitude  data  on  an  airfoil  pitching  up  to  60°.  Clearly,  an  extensive  amount  of 
experimental  effort  needs  to  be  performed  to  obtain  a  clear  understanding  of  the  fluid 
dynamic  mechanisms  precipitating  the  occurrence  of  the  dynamic  vortices  associated 
with  the  unsteady  flowfields  surrounding  airfoils  pitching  at  high  rates  to  large 
amplitudes.  This  experimental  effort  intends  to  add  to  that  understanding  by  exploring 
a  well  defined  but  limited  parameter  range  of  pitch  rates  and  Reynolds  numbers. 

II.  Experimental  Arrangement 

All  data  were  obtained  in  the  USAF  Academy  Aeronautics  Laboratory  low-speed, 
2  ft  X  3  ft  subsonic  wind  tunnel.  The  model  used  was  a  full  span,  6  inch  chord,  NACA  001 5 
airfoil.  The  airfoil  was  instrumented  with  eighteen  Endevco  8507-2  miniature  pressure 
transducers  close-coupled  to  the  surface  ports,  as  shown  in  Fig.  1.  The  signals  from  each 
transducer  were  amplified  through  Dynamics  7512B  amplifiers.  The  pitching  motions  of 
the  airfoil  were  generated  by  a  PDP  11/03  computer,  which  commanded  a  Control 
Systems  Research  Index-Syn  stepper  motor/driver  system.  This  is  a  modified  version  of  the 
two-degree-of-freedom  oscillator  developed  by  Francis,  et.  al.  (Ref.  16).  The  PDP  1 1/03 
was  slaved  to  a  PDP  1 1/45,  which  performed  the  supervisory  control  and  data  acquisition. 
Since  the  airfoil  motions  were  very  repeatable,  all  eighteen  pressure  ports  were  located 
on  one  side  of  the  airfoil  for  better  data  resolution.  This  resulted  in  some  minor  errors  in 
the  data  between  the  upper  and  lower  surfaces  due  to  difficulties  in  returning  the  tunnel 
to  the  precise  flow  speed  for  the  opposite  surface  after  turning  the  airfoil  over.  The 
reversal  of  the  airfoil  was  necessitated  by  the  fact  that  the  stepping  motor  would  turn  in 
only  one  direction. 
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Figure  1.  Pressure  Port  Locations 


Experiments  were  conducted  employing  a  variety  of  constant  rate  pitching 
motions  and  flow  speeds.  Table  1  shows  the  experimental  matrix  with  a  given  in  deg/sec 
and  L/a  in  ft/sec,  and  with  the  interior  values  given  by  the  non-dimensional  pitch  rate 
parameter  a  +  =  a  c/L/«.  All  of  the  motions  were  performed  by  pitching  the  airfoil  from 
angles  of  attack  of  zero  degrees  to  sixty  degrees  and  then  stopping  the  motion.  See 
Fig.  2. 


Table  1 
a (o/s) 


U.(f/s) 

230 

460 

920 

1380 

20 

0.1 

0.2 

0.4 

0.6 

40 

— 

0.1 

0.2 

— 

60 

— 

— 

— 

0.2 

80 

— 

0.05 

0.1 

— 

The  Reynolds  number  based  on  chord  varied  from  47500  to  190000  The  major 
objectives  of  these  experiments  were  to  examine  the  effects  of  varying  pitch  rate  at 
constant  Reynolds  number,  varying  Reynolds  number  at  constant  pitch  rate,  and  holding 
a-*-  constant  while  varying  both  pitch  rate  and  Reynolds  number  For  the  nine  cases 
shown  in  Table  I,  data  were  taken  for  both  surfaces  of  the  airfoil  by  performing  twenty- 
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Figure  2.  Angle  of  attack  versus  time 


five  repetitions  of  the  pitching  motions,  taking  200  data  points  per  sensor,  and  taking 
the  ensemble  average. 

III.  Experimental  Results 

The  objectives  of  these  experiments  were  to  examine  the  unsteady  aerodynamics 
surrounding  a  pitching  airfoil  and  the  time  dependent  fluid  mechanics  associated  with 
varying  pitch  rate,  Reynolds  number,  and  the  non-dimensional  pitch  rate  parameter,  a+. 
The  matrix,  shown  in  Table  1  above,  was  designed  to  accomplish  this  task  within  the 
experimental  constraints  of  the  control,  instrumentation,  and  data  acquisition  system. 
The  pressure  data  were  used  to  calculate  pressure  coefficients  which  were  in  turn 
integrated  to  compute  lift  coefficients.  Both  types  of  curves  were  plotted  in  the  usual 
way;  Cp  vs  x/c  and  C/  vs  a.  Since  these  experiments  were  limited  to  a  single  airfoil, 
varying  the  Reynolds  number  was  synonymous  with  varying  the  flow  velocity  in  the 
tunnel. 

A.  Section  Lift  Coefficient  versus  Alpha 

The  effect  of  varying  the  pitch  rate  while  holding  Reynolds  number  constant  was 
dramatic,  particularly  at  the  lower  flow  velocities.  At  Re  =  47500,  was  a  strong 
function  of  pitch  rate,  varying  from  2  to  4.3  as  pitch  rate  was  increased  from  230  deg/s  to 
1380  deg/s.  At  this  Reynolds  number,  dCi/da  was  also  a  strong  function  of  pitch  rate. 
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varying  from  4.3/rad  to  8/rad.  As  Reynolds  number  increased  to  95000  and  190000, 
however,  their  dependence  is  less  pronounced.  Fig.  3  shows  C;  vs  a  for  various  pitch  rates 
at  three  different  Reynolds  numbers,  and  Fig  4  shows  the  effect  on  Cimax  oi  increasing 
pitch  rate  at  constant  Re. 


a  (d«g] 


a  (dag) 


0  10  20  30  40  30  60 

a  (deg) 

Figure  3.  Ci  vs  a  for  constant  Uxat  various  pitch  rates 


Fig  4  shows  (from  Figure  3(a)  only)  a  somewhat  linear  increase  m  with  respect  to 


pitch  rate.  Gormont  (Ref.  17)  and  Strickland  and  Graham  (Ref.  18)  have  shown  a  square 


root  dependency  on  non-dimensional  pitch  rate  for  inception  of  dynamic  stall.  Gormont 
also  shows  a  similar  correlation  for  lift  stall  Also,  Walker,  et.  al.  (Ref  14)  and  Helin  and 
Walker  (Ref  15)  presented  near-surface  hot-wire  data  that,  while  being  impossible  to 
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Figure  4.  Cimax  vs  a  at  various  Reynolds  number 

calculate  lift  coefficients  from,  indicates  support  for  this  conclusion.  While  for  steady 
flows,  a  direct  correlation  exists  between  C/max  arid  stall,  interestingly,  for  these  types  of 
unsteady  flows,  the  correlation  is  not  straightforward.  In  addition  to  the  effect  of  pitch 
rate  on  C/max.  there  is  a  general  dependence  of  the  slope  of  the  lift  curve  and  of  the 
angle  of  attack  at  which  C/^ax  occurs;  and  again,  this  dependence  is  less  pronounced  as 
Reynolds  number  increases. 

The  effect  of  varying  Reynolds  number  while  holding  pitch  rate  constant  is  rather 
the  opposite.  C/m^x  reduced  significantly  when  increasing  Reynolds  number  from 
47500  to  95000  or  1425000,  but  is  only  slightly  affected  by  increasing  Re  from  95000  to 
190000.  Fig.  5  shows  C/  vs  a  for  various  Reynolds  numbers  at  three  different  pitch  rates, 
and  Fig.  6  shows  the  effects  of  C/max  of  increasing  Re  at  constant  pitch  rate.  The  angle  of 
attack  at  which  C/^ax  occurs  is  also  inversely  dependent  on  Reynolds  number,  as  is  dQ/da 
for  rotation  rates  of  920  an  '  1380  deg/s.  Oddly  enough,  however,  the  lift  curve  slopes  at 
the  pitch  rate  of  460  degis  are  all  much  the  same. 
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The  effect  of  maintaining  the  non-dimensional  pitch  rate  parameter,  a+,  constant 
within  the  range  of  Reynolds  numbers  investigated  here  is  of  special  significance.  Two 
cases  were  studied:  that  of  a  +  =  0.1  and  q  +  =  0.2.  As  can  be  seen  from  Table  1,  these 
were  obtained  by  increasing  both  pitch  rate  and  flow  velocity  proportionately.  Figs.  7 
and  8  show  C/  vs  a  curves  for  the  respective  constant  cases. 


Figure?.  a+  =0.1, a  =  230, 460, 920 deg/sec,  U«  =  20, 40, 80 ft/sec 


Figure  8.  a  ♦  =  0.2,  a  =  460, 920, 1 380  deg/sec,  U«  =  20, 40,  80  ft/sec 
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As  one  can  readily  see,  each  set  of  constant  a  +  curves  are  very  similar.  For  a+  =  O.t, 
each  of  the  three  curves  start  out  at  a  high  initial  slope  for  the  first  two  or  three  degrees, 
then  proceed  linearly  at  a  slope  of  approximately  4.2  until  reaching  a  Qmax  about  2.0 
at  31  or  32  degrees.  The  higher  pitch  rates  produce  slightly  higher  values  of 
slightly  higher  angles  of  attack,  but  to  this  point  they  are  virtually  the  same  within  the 
error  band  of  the  data.  The  curves  then  drop  off  rapidly  -  the  higher  the  rotation  rate, 
the  sharper  the  drop  with  respect  to  angle  of  attack  -  and  oscillate  about  a  C/  of  1.2  out  to 
60 degrees.  The  three  lift  curves  at  a*  =  0.2  do  not  behave  in  quite  the  same  manner  as 
those  at  a  +  =  0.1.  They  start  out  for  the  first  20  degrees  or  so  with  the  same  slopes  as 
one  another  but  not,  as  with  the  lower  value,  linear.  They  then  deviate  somewhat  from 
one  another,  peaking  out  at  the  same  value  of  C/max  of  2.2,  but  at  angles  of  attack 
varying  from  31  to  41  degrees.  There  seems  to  be  no  specific  trend  which  can  be 
attributed  to  either  rotation  rate  or  flow  velocity,  since  the  curve  that  appears  to  be  in 
the  middle  is  actually  that  for  the  lowest  pitch  rate.  At  these  high  pitch  rates  this  minor 
anomaly  could  be  attributed  to  experimental  error.  The  conclusion  remains,  however, 
that  the  non-dimensional  pitch  rate  parameter,  a+.  is  at  this  point,  the  major 
determining  factor  for  the  shape  and  magnitude  of  the  lift  curve.  Thus,  since  curves  of 
constant  a+  are  nearly  similar,  examing  the  effects  of  changing  a+  is  closely  akin  to 
analyzing  the  effects  of  changing  the  pitch  rate  at  constant  Reynolds  number  as  in  Fig.  3. 
One  can  see  from  there  and  from  Figs.  7  and  8  that  increasing  a+  results  in  increases  in 
dCiida,  Cimax  .and  a  at  Cimax  In  addition,  lower  values  of  a  +  result  in  more  abrupt 
decreases  in  Q  occurring  more  rapidly  after  having  reached  Q^ax 
B.  Pressure  Coefficient  vs  x/c 

The  Appendix  contains  Figs.  A.  1  (a),  (b),  (c),  (d)  through  A. 9  (a),  (b),  (c),  (d).  The 
numbers  1-9  correspond  to  each  of  the  cases  shown  in  Table  1  starting  with 
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C/oo  =  20f\Jsec,  d  =  230deglsec,  a+  =0.1,  following  left  to  right  and  ending  with 
Ua>  =  SOWsec,  a  =  920deglsec,  a+ =  0.1.  The  (a)  -  (d)  discriminants  refer  to  angles  of 
attack  of  15°,  30°,  45°,  and  60°  respectively. 

Walker  et.  al.  (Ref.  14)  presented  qualitative  data  using  a  smoke-wire  flow 
visualization  technique  and  near-surface  single  element  hot-film  velocity  magnitude 
measurements  for  two  cases  essentially  the  same  as  two  of  the  cases  examined  herein; 
a  +  =0.2  and  0.6.  In  these  experiments  they  observed  strong  vortices  forming  above  the 
suction  surface  of  the  airfoil  due  to  the  shear  layer  interaction  between  the  freestream 
and  the  leading  edge  separation  bubble,  formed  as  the  airfoil  pitches  up  and  to  the 
added  rotational  energy  imparted  to  the  flow.  They  were  able  to  correlate  the 
movements  of  these  vortices  with  velocity  magnitude  peaks  measured  with  the  hot-film 
sensors.  Figs.  A.2  and  A.4  show  the  corresponding  surface  pressure  measurement 
experiments  of  these  two  cases;  the  upper  surface  pressure  peaks  are  in  good  agreement 
with  the  spatial  positions  of  the  vortices. 

Fig.  A.1  shows  the  lowest  pitch  rate  examined  corresponding  to  an  a+  of  0.1. 
Except  for  the  fact  that  the  NACA  0015  airfoil  would  be  stalled  at  15  degrees  in  steady 
flow,  (a)  shows  a  rather  ordinary  Cp  plot  for  high  angles  of  attack,  with  the  exception  of 
the  first  20%  of  the  upper  surface.  Here  one  can  observe  a  pressure  peak  at  the  16% 
chord  point  indicating  the  presence  of  the  leading  edge  vortex.  Also,  there  is  no  high 
suction  peak  very  near  the  leading  edge  that  is  present  in  steady  flow.  This  is  true  for  all 
of  the  other  cases  studied  with  the  exception  of  the  lowest  (  a+ =  0.05)  case.  This  could, 
in  part,  be  due  to  the  first  sensor  location  but  is  more  likely  due  to  the  interaction  of  the 
leading  edge  vortex  reverse  flow  in  the  vicinity  of  the  leading  edge  and  the  freestream 
flow  producing  locally  reduced  flow  velocities.  In  (b)  we  can  see  that  the  vortex  is  still 
producing  a  great  deal  of  lift  and  has  traversed  to  a  point  slightly  aft  of  mid-chord  In  (c) 
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we  can  now  tell  that  the  vortex  has  separated  from  the  airfoil  causing  it  to  stall. 

Increasing  a+  to  a  value  of  0.2,  Fig.  A.2  (a)  shows  characteristics  similar  to  A.1  (a) 
but  higher  in  value;  and  then  (b)  shows  a  very  high  pressure  peak  of  almost  -6.0  at  24% 
chord  indicating  the  presence  of  a  stronger  vortex  than  previously  observed.  By  the  time 
the  airfoil  reaches  45°  in  (c),  the  vortex  has  moved  past  mid-chord  and  off  the  airfoil, 
substantially  reducing  lift.  At  60°,  the  airfoil  is  completely  stalled.  Also,  an  interesting 
pressure  anomaly  has  entered  the  picture. 

From  8%  to  56%  chord  at  an  angle  of  attack  of  45°,  the  lower  surface  pressure 
coefficient  is  slightly  greater  than  one.  This  could  be  attributed  to  experimental  error 
except  for  two  elements:  1)  Francis,  et.  al.  (Ref.  1 1)  noticed  this  sort  of  phenomenon  in 
their  experiments  and  stated  that  it  should  be  investigated  further;  2)  in  these  sets  of 
experiments,  this  phenomenon  is  observed  to  become  more  pronounced  at  higher  angles 
of  attack  and  at  higher  values  of  a+  .  It  appears  to  occur  at  q+  levels  of  0.2  or  greater, 
where  a  great  deal  of  rotational  energy  relative  to  the  freestream  is  injected  into  the 
flow  producing  large  flow  curvatures.  King  (Ref.  19)  has  pointed  out  that  a  Wr  term 
from  Crocco's  theorem  as  presented  in  Liepmann  and  Roshkow  (Ref.  20)  may  be  the 
culprit  by  increasing  the  local  stagnation  pressure. 

For  the  case  of  q-*-  =  0.4,  Fig.  A. 3  (a)  indicates  a  very  high  value  of  Cp  over  the 
entire  airfoil  similar  in  nature  for  the  previous  case  but  greater  in  magnitude.  At  an 

1 

angle  of  attack  of  30°,  (b)  shows  a  very  targe  peak  in  Cp  of  greater  than  -9.0  at  the  10%  j 

( 

chord  point.  By  the  time  the  airfoil  has  reached  an  angle  of  attack  of  45°  in  (c),  the  vortex 
has  moved  downstream  to  the  34%  chord  point  and  away  from  the  airfoil.  At  60°,  shown 
in  (d),  the  vortex  has  moved  further  downstream  but  is  still  exhibiting  a  great  deal  of  ( 

j 

influence  on  the  flowfield  around  the  airfoil.  l 

{ 

i 

Fig.  A.4  depicts  pressure  coefficients  for  the  last  of  the  four  cases  at  Um  =  20  ftlsec.  I 

i 
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This  is  the  case  of  the  highest  a+  examined.  At  a+  =  0.6  the  airfoil  is  pitching  up  at  a 
rate  of  1380  deg/sec.  At  15°  in  (a),  Cp  curves  look  about  the  same  as  before  except  a  bit 
higher.  In  (b),  the  vortex  hasn't  really  formed  yet,  but  the  pressure  peak  at  6%  chord  is  at 
-9.8  and  the  flow  is  still  attached  over  the  entire  airfoil.  In  addition,  the  maximum  lower 
surface  pressure  coefficient  has  reached  1.5.  By  45°  the  vortex  has  fully  formed, 
producing  a  pressure  peak  greater  than  -14.0  at  the  14%  chord  point.  These  conditions 
have  also  resulted  in  a  maximum  lower  surface  Cp  of  2.0  at  47%  chord.  As  the  airfoil 
reaches  60“  in  (d),  the  vortex  has  moved  to  24%  chord  As  the  airfoil  reaches  60“  in  (d), 
the  vortex  has  moved  to  24%  chord  and  away  from  the  surface  resulting  in  a  reduced  Cp 
of -9.0. 

Figs  A.5  and  A.6  show  pressure  coefficients  for  two  pitch  rates  at  Ua,=  40  ft/sec 
corresponding  to  non-dimensional  pitch  rates  of  0.1  and  0.2  respectively.  They  are 
remarkably  similar  to  Figs.  A.1  and  A. 2  with  only  minor  differences  in  peak  values  and 
locations. 

The  only  data  taken  at  60  ft/sec  are  for  a  +  =  0.2  with  a  corresponding  pitch  rate  of 
1380  deg/sec  shown  in  Fig.  A.7.  Again,  as  with  Fig.  A.6,  the  pressure  coefficient  data  is 
very  similar  to  that  of  Fig.  A.2  -  the  only  exception  being  that  of  A.7  (c),  where  a  relatively 
high  peak  is  seen  at  34%  chord. 

Two  final  sets  of  experimental  data  were  taken  at  Um  =  80  ft/sec.  These  are  shown 
in  Figs.  A. 8  and  A.9  for  a  +  values  of  0.05  and  0. 1  respectively.  In  A.8  (a)  we  see  Cp  curves 
much  like  those  of  unstalled  high  angles  of  attack  for  this  airfoil.  At  30“  it  is  apparent 
that  the  airfoil  has  stalled  completely,  indicating  that  the  leading  edge  vortex  has  either 
not  formed  or  does  not  have  the  requisite  strength  to  produce  the  strong  reverse  flow 
velocities  observed  previously. 

Figs.  A.9  (a),  (c),  and  (d)  look  very  much  like  their  counterparts  in  Figs.  A.2  and  A, 5. 
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A.9  (b),  however,  is  somewhat  different.  Whereas  the  first  two  cases  where  a  +  =  0.1  at 
anangleof  attack  of  30°  produced  suction  surface  pressure  coefficients  of  -1.5  to -1.8 
near  the  leading  edge,  increasing  to  -3.5  at  47%  to  53%  chord,  and  dropping  to  -0.5  at 
85%  chord,  this  one  starts  at  -3.0  near  the  leading  edge,  peaks  at  -4.0  at  34%  chord,  and 
falls  rapidly  to  -0.5  at  65%  chord.  While  this  will  produce  almost  the  same  lift 
coefficients,  the  indication  in  this  case  is  that  the  leading  edge  vortex  develops 
somewhat  more  rapidly  and  is  somewhat  stronger  than  in  the  previous  a+  =0.1  cases. 
Additionally,  a  greater  portion  of  the  aft  section  of  the  airfoil  is  separated. 

IV.  Conclusions 

The  experimental  pressure  measurements  presented  herein  have  shown 
relationships  between  the  unsteady  lifting  forces  on  a  pitching  airfoil  and  pitch  rate, 
Reynolds  number,  and  non-dimensional  pitch  rate,  a+  .  The  current  investigations  were 
limited  to  the  effects  of  these  parameters  on  pressure  and  lift  coefficients.  With  regard 
to  lift  coefficients  it  was  seen  that,  for  constant  Reynolds  number  (or  flow  velocity), 
increasing  pitch  rate  (or  non-dimensional  pitch  rate)  increased  C/max.  dC/da  up  to  C/max. 
and  in  general,  the  angle  of  attack  which  C/max  occurs.  These  effects  were  more  dramatic 
at  lower  Reynolds  numbers.  Increasing  Reynolds  number  while  holding  pitch  rate 
constant  had  the  opposite  effect.  The  effect  of  holding  the  non-dimensional  pitch  rate 
constant  while  allowing  pitch  rate  and  Reynolds  number  (flow  velocity)  to  vary  is  of 
special  interest.  For  both  cases  studied,  a+  =  0. f  and  0.2,  the  C/ vs  a  curves  were  very 
similar  at  constant  a  +  -  only  the  angle  of  attack  for  C/^ax  showed  inconsistencies  at  the 
higher  value. 

Examination  of  the  non-integrated  parameters  or  pressure  coefficients,  proved 
more  enlightening.  Not  only  did  the  Cp  vs  x/c  curves  show  the  same  overall  effects  as  the 
C/  vs  a  curves  did,  but  they  showed  the  same  effects  at  discrete  points  on  the  airfoil. 
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Pressure  coefficients  at  given  angles  of  attack  for  constant  a+  values  showed  the  same 
correlation  as  the  lift  coefficients.  In  addition,  the  movement  of  the  "dynamic  stall" 
vortices  and  their  effects  on  the  airfoil  under  different  conditions  could  be  determined. 
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Undergraduate  Airbreathing  Propulsion  At  The 
US  Air  Force  Academy 

Jack  D.  Mattingly*  and  William  H.  Heiser** 

Abstract 

This  paper  summarizes  the  status  of  the  teaching  of  airbreathing  propulsion  to  the 
undergraduate  students  at  the  U.S.  Air  Force  Academy.  This  is  one  of  the  most  active  and 
productive  propulsion  curricula  in  the  country.  The  emphasis  of  this  article  is  on  course 
philosophy  and  content  in  general,  and  on  the  goals  and  experiences  of  the  flagship 
course  on  propulsion  design  in  particular. 

I.  Introduction 

A  very  active  curriculum  in  airbreathing  propulsion  exists  in  the  Department  of 
Aeronautics  of  the  U.S.  Air  Force  Academy,  for  a  number  of  mutually  supportive  reasons. 
To  begin  with,  many  cadets  attend  the  Academy  in  order  to  become  pilots;  therefore, 
they  have  a  natural  curiosity  about  and  a  vital  interest  in  how  jet  engines  work. 
Furthermore,  the  cadets  realize  that  such  knowledge  will  be  necessary  because  their 
careers  are  likely  to  include  one  or  more  interesting  and  important  assignments  in 
laboratories,  procurement  offices,  maintenance,  ground  and  flight  test,  or  as  staffs  of 
the  Commands  and  Pentagon.  Finally,  unusually  heavy  support  for  the  curriculum  is 
available  in  the  form  of  a  faculty  rich  in  experience  in  the  development  and  operation  of 
aircraft  engines,  generous  supplies  of  computer  time  and  codes,  and  immediate  access  to 
excellent  laboratory  facilities. 

In  fact,  the  Aeronautics  Laboratory  may  be  the  only  place  in  the  country  with  enough 
financial  and  logistic  support  to  allow  undergraduate  students  to  routinely  experiment 
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with  real  jet  engines.  The  laboratory  has  four  engine  test  cells  (see  Table  1).  Three  are 
dedicated  to  the  testing  of  jet  engines  and  one  is  dedicated  to  solid  propellant  rocket 
thrust  tests  and  internal  combustion  engine  testing. 


Table  1  Aeronautics  Laboratory  Propulsion  Equipment 


Test  Cell  # 

Description 

1 

One  J85  turbojet  engine  -2800  lb  thrust,  non-  augmented 

2,3 

Two  J69  turbojet  engines -1000  lb  thrust  (one  is  planned  to 
be  replaced  by  an  F109  turbofan  engine  - 1300  lb  thrust) 

4 

5'  diameter  solid  propellant  rocket  motor  (500  lb  thrust, 

3  sec  burn  time,  C*  =  5000  ft/sec) 

All  cadets  graduating  from  the  Academy  (approximately  1,000/year)  perform  a  static 
thrust  test  with  the  J69  turbojet  engines  as  part  of  the  'Engineering  Thermodynamics" 
course  (a  required  course  for  all  graduates)  taught  by  the  Department  of  Aeronautics. 
Since  the  students  are  interested  in  and  curious  about  the  operation  of  airbreathing 
engines  and  rocket  nozzles,  this  course  demonstrates  the  application  of  engineering 
thermodynamics  by  analyzing  these  machines.  The  static  thrust  test  provides  the 
opportunity  to  compare  analytical  and  experimental  results. 

The  result  of  these  forces  is  a  successful  aeronautical  engineering  program  which 
educates  approximately  80  new  cadets  each  year,  and  was  the  breeding  ground  (in  1976) 
for  the  famous  "Oates  Notes"  (Ref.  1),  which  have  become  the  bestselling  AIAA 
Education  Series  Book,  The  Aerothermodvnamics  of  Gas  Turbine  and  Rocket  Propulsion  . 
by  Gordon  C.  Oates. 

II .  The  Curriculum 

The  aeronautical  engineering  major  is  currently  composed  of  thirty-three  semester 
hours  of  major  courses  in  addition  to  the  1 1 1  semester  hours  contained  in  the  Academy's 
Core  Curriculum.  Because  of  the  limited  number  of  courses  in  the  major,  topics  that 
typically  would  be  taught  in  separate  courses  are  combined  and  split  between  courses  to 
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provide  the  required  coverage.  For  example,  gas  dynamics  is  incorporated  into  both  the 
Propulsion  I  course  (one-dimensional  flow)  and  the  Aerodynamics  II  course  (two- 
dimensional  flow).  The  curriculum  (Table  2)  provides  the  opportunity  for  student  choice 
of  design  course  (  6  semester  hour  course)  and  one  Aero  option.  Thus  the  student  can 
emphasize  a  specific  area  of  interest,  such  as  flight  mechanics,  propulsion,  aerodynamics, 
or  structures.  j 

_ Tabic  2  Aeronautical  Engineering  Major _ 

•  1 1 1  Semester  hours  of  core  courses 
’  33  Semester  hours  of  major  courses 

Junior  Year 

•  Flight  Mechanics  I 

-  Propulsion  I 

•  Aerodynamics  I 

•  Aerodynamic  Structures 

•  Analytical  &  Numerical  Methods 
in  Aeronautics 

•  Senior  Year 

•  Aeronautical  Laboratory 

•  Aerothermodynamics 

•  Aerodynamics  II 

•  Design  Option  •  Aircraft  Design  or 

Propulsion  Design 

•  Aero  Option  •  Flight  Mechanics  II 

Propulsion  II 
Flight  Test  Techniques 
Special  Topics 
Independent  Study 


The  "propulsion  track"  is  the  equivalent  of  an  undergraduate  minor  in  airbreathing 
propulsion.  It  consists  of  three  courses,  usually  taken  in  sequence,  which  carry  the 
student  from  the  fundamental  principles  through  to  the  complete  design  of  an  engine 
which  can  accomplish  a  very  specific  mission.  The  three  courses  will  now  be  described  in 
their  natural  order. 
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1.  Propulsion  I 

This  introductory  course  is  required  in  order  to  obtain  an  Aeronautics  degree  at 
*^he  Academy.  It  equates  to  a  three  unit  course  in  the  semester  system,  and  includes  the 
aforementioned  80  students  annually.  The  emphasis  of  the  course  is  upon  the 
aerodynamic  and  thermodynamic  behavior  of  the  one-dimensional  flow  of  calorically 
perfect  gases.  The  central  feature  of  the  course  is  the  development  of  methods  to  relate 
the  overall  performance  to  internal  design  parameters  of  rocket  and  gas  turbine  engines, 
in  the  manner  of  Oates  (Ref.  1).  This  feature  receives  special  emphasis  because  every 
student  participates  in  one  rocket  firing  and  one  gas  turbine  engine  run  where 
experimental  data  is  systematically  compared  with  predictions.  The  influence  of  the 
propulsion  system  and  the  proper  engine  cycle  on  aircraft  performance  are  learned 
through  a  design  project.  For  the  repetitive  calculations  required  in  the  design  project, 
students  use  computer  codes  based  on  the  on-design  and  off-design  cycle  analysis 
presented  in  class. 

Further  insight  into  the  content  and  flow  of  this  course  can  be  obtained  by  looking  at 
the  formal  table  of  contents  below. 

One  Dimensional  Flows 

Review  of  fundamental  laws 
Normal  and  oblique  shock  waves 
Duct  flow  with  heating,  friction,  and  area  change 
Rocket  Propulsion 

Performance  definition  and  prediction 
Nozzle  theory  and  design 
Rocket  motor  laboratory  test 
Aircraft  Propulsion 

General  concepts  (energy  chain  and  installed  vs.  uninstalled  thrust) 

Design  point  cycle  analysis; 
ideal  ramjet 
ideal  turbojet 

ideal  afterburning  turbojet 
ideal  turbofan 
real  turbojet 

real  afterburning  turbojet 
real  turbofan 
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sensitivity  analysis 
optimum  performance 
Off-design  cycle  analysis: 
real  turbojet 
real  turbofan 
Components: 
inlets 

compressors 

burners 

turbines 

afterburners 

nozzles 

Turbine  engine  laboratory  test 
Cycle  design  project 

2.  Propulsion  II 

This  course  also  equates  to  three  units  in  the  semester  system,  and  includes  about 
20  students  annually.  The  emphasis  is  upon  qualitative  and  quantitative  aspects  of  the 
design  of  the  major  components  of  turbine  engines.  For  example,  methods  are 
developed  which  enable  the  design  of  compressor  and  turbine  stages  and  inlet 
configurations,  as  well  as  the  calculation  of  the  products  of  combustion  reactions.  Also 
included  is  the  analysis  of  more  advanced  jet  engine  cycles,  and  a  final  project  which  tests 
the  understanding  of  several  of  the  course  topics  at  once. 

Although  this  is  generally  considered  as  preparation  for  the  design  course  which 
follows,  it  need  not  necessarily  be  taken  for  that  purpose  Some  students  terminate  with 
this  course,  a  few  others  take  it  simultaneously  with  the  design  course,  and  others 
occasionally  go  on  to  the  design  course  without  it.  Further  appreciation  of  the  scope 
should  be  made  possible  by  the  formal  table  of  contents  below. 

Compressor  Design 

Cascade  aerodynamics 
Stage  selection 
Radial  equilibrium 
Multi-staging 
Reduced  compressor  map 
Post-stall  behavior 

Turbine  Design 

Cascade  aerodynamics 
Work  relationships 
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Structural  analysis 

•  Component  Matching  and  Sizing 

•  Combustion  Chemistry 

Mixtures  of  perfect  gases 

Concepts  (heat  of  reaction,  equilibrium,  law  of  mass  action) 
Chamber  condition  calculations 
Equilibrium  and  frozen  nozzle  flows 
Chemical  kinetics 

•  Stationary  Components 

Inlets 

Burners 

Mixers 

Afterburners 

Nozzles 

•  Advanced  Cycles 

Turbofan  with  mixing  and  afterburning  (on-  and  off-design) 
Thrust  augmentation  devices 
Component  Design  Project 


3.  Propulsion  Design 

This  is  both  the  centerpiece  and  the  culmination  of  the  "propulsion  track."  It  is  a 
double-unit  course,  equaling  two  three-unit  courses  in  the  semester  system,  and 
including  about  20  students  annually.  In  contrast  to  the  two  preceding  courses,  which 
strongly  emphasize  basic  concepts  and  specific  analytical  tools,  this  course  contains  many 
features  which  are  intended  to  reproduce  the  design  experience  via  the  design  of  an 
engine. 

There  is  no  absolute  roadmap  for  the  design  of  a  gas  turbine  engine.  The  steps 
involved  depend,  for  example,  on  the  experience  of  the  company  and  the  people 
involved,  as  well  as  on  the  nature  of  the  project.  A  revolutionary  new  engine  will  require 
more  analysis  and  iteration  than  the  modification  of  an  existing  powerplant. 
Nevertheless,  there  are  generalized  representations  of  the  design  process  which  are 
informative  and  useful.  One  of  these,  which  depicts  the  entire  development  process,  is 
shown  in  Fig.  1  which  is  taken  from  Gas  Turbine  Theory  by  Cohen,  et  al  (Ref.  2).  The 
portion  of  Fig.  1  enclosed  by  the  dashed  line  is  of  paramount  importance  here  because 
that  is  the  territory  covered  by  this  course.  Although  the  boundary  can  be  somewhat 
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Figure  1.  Gas  Turbine  Engine  Design  System 

altered,  for  example  by  including  control  system  studies,  it  can  never  encompass  the 
hardware  phases,  such  as  manufacturing  and  testing. 

We  would  like  to  make  clear  what  we  believe  is  special  about  the  design  process,  for 
that  is  what  this  course  attempts  to  make  happen.  Every  designer  has  a  mental 
description  of  what  designing  really  is,  so  our  version  is  not  likely  to  be  exhaustive. 
Nevertheless,  the  list  below  contains  critical  elements  with  which  few  would  disagree. 


The  design  process  is  both  started  by  and  constrained  by  an  identified  need. 
Although  many  legitimate  solutions  exist,  none  can  be  identified  as  unique  or 
optimum.  The  final  selection  always  involves  judgment  and  compromise. 

In  the  case  of  systems  design,  such  as  aircraft  and  engines,  too  many  legitimate 
solutions  exist.  Systematic  methods  must  be  found  to  home  in  on  the  best. 
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The  process  is  inherently  iterative,  often  requiring  the  return  to  an  earlier  step 
when  prior  assumptions  lead  to  a  dead  end. 

Many  technical  specialties  are  interwoven.  For  example,  gas  turbine  engine 
design  involves  at  least  thermodynamics,  aerodynamics,  heat  transfer, 
combustion,  structures,  materials,  manufacturing  processes,  instrumentation, 
and  controls. 

Above  all,  the  design  of  a  complex  system  requires  active  participation  and 
disciplined  communication  by  every  one  involved.  Because  each  part  of  the 
system  influences  all  the  others,  the  best  solutions  can  only  be  discovered  (and 
major  problems  and  conflicts  avoided)  if  the  participants  share  their  findings 
clearly  and  regularly. 

Some  of  the  special  features  we  have  found  that  greatly  contribute  to  the 
accomplishment  of  these  goals  include: 

■  Basing  the  design  on  a  contemporary  Request  for  Proposal  (RFP)  provided  by  an 
interested,  outside  "sponsor",  such  as  the  Air  Force  Flight  Dynamics  Laboratory 
(AFFDL)  or  the  Air  Force  Air  Staff.  Figure  2  shows  the  mission  of  a  recent  project 
sponsored  by  AFFDL. 

Working  in  concert  with  an  aircraft  design  course  on  the  same  RFP,  in  which 
case  they  play  the  role  of  an  airplane  company.  Frequent  contact  and 
negotiation  between  courses  helps  improve  results,  while  providing  a  lot  of 
reality. 

Supplying  very  user  friendly  computer  codes  for  engine  on-  and  off-design  cycle 
analysis,  along  with  liberal  amounts  of  computer  time,  so  that  the  students  can 
concentrate  on  the  real  objective  of  finding  the  "best"  solutions. 
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Figure  2.  Air-to-Air  Fighter  Mission  Profile 

•  Requiring  individual  design  of  specific  components  (e.g.  fan,  burner,  turbine, 
nozzle)  along  with  correct  aerothermodynamic,  mechanical  and  geometrical 
inter-relationships.  This  emphasizes  the  need  for  frequent,  disciplined 
communications,  and  culminates  in  a  single  design  drawing  of  the  integrated 
engine  as  a  proof  test. 

•  Recycling  results  when  necessary  to  reach  the  "best"  solution. 

Scheduling  periodic  written  and/or  verbal  reports  in  order  to  insure  the  flow  of 
information,  and  to  sharpen  communication  skills. 

•  Rotating  leadership  roles  (e  g.  program  manager,  chief  engineer,  vice  president 
of  engineering)  among  students  in  order  to  familiar'ie  them  with  the 
experiences  of  organizing  and  managing. 
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The  design  process  as  modified  for  this  course  is  shown  in  Fig.  3  and  the  formal  course 
content  which  is  intended  to  achieve  these  goals  is  as  follows; 

•  Engine/Aircraft  System  Constraint  Analysis 

Takeoff  thrust-to-weight  (thrust  loading) 

Takeoff  weight-to-wing  area  (wing  loading) 

•  Engine/Aircraft  System  Mission  Analysis 

Takeoff  thrust 
Takeoff  weight 
Takeoff  wing  area 
Specific  excess  power  (Pj) 

•  Engine  Design  Point  Analysis 

Ranges  of  design  parameters 
Initial  performance  estimates 
Engine  Off-Design  Analysis 

Selection  of  design  parameters 
Installation  penalties 
Sizing 
Throttling 

Final  performance  analysis 
Takeoff  thrust 
Takeoff  weight 
Internal  flow  properties 

•  Component  Design 

Aerothermodynamic 

Structural 

•  Subsystem  Design 

•  Engine  Layout 

•  Reports  and  Presentations 


We  have  come  to  recognize  two  events  that  are  somehow  intimately  connected  to 
the  success  of  Propulsion  Design.  First,  the  emotional  rewards  are  highest  when  it 
appears  that  a  successful  engine  cannot  be  found  at  some  point  during  the  design  cycle. 
When  the  barriers  (not  the  rules)  are  legitimately  broken  by  the  usual  method  of 
exhaustively  searching  for  innovative  solutions,  the  loose  coalition  of  frustrated  students 
is  transformed  into  a  confident  team.  Moreover,  the  lesson  has  real  applicability  to  their 
prospective  careers  in  engineering.  Second,  there  seems  always  to  be  at  least  one 
moment  at  which  the  entire  class  makes  a  mistake  of  unimaginably  large  proportions. 
Examples  of  this  have  included  area  mismatches  between  components,  designing  an 
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engine  with  the  airflow  intended  for  two,  or  the  compressor  and  turbine  attached  to  the 
same  shaft  having  different  rotational  speeds.  These  episodes  have  a  very  sobering 
effect  on  the  students  and,  after  the  recriminations  dwindle  away,  result  in  a  far  greater 
appreciation  of  the  vital  need  for  frequent,  accurate,  and  intelligent  communication. 

Through  the  entire  design  process,  it  is  less  important  that  specific  milestones  be 
reached  according  to  a  preset  schedule  than  that  the  students  achieve  comprehension  of 
each  successive  step,  and  that  they  believe  in  their  ownership  of  the  design  decisions. 
This,  above  all  things,  makes  their  experience  take  life. 


The  propulsion  curriculum  has  been  based  on  the  analytical  methods  for  cycle  and 
component  analysis  as  set  down  by  Gordon  C.  Oates  in  "Oates  Notes"  during  1975-76 
while  a  Distinguished  Visiting  Professor  with  the  Department  of  Aeronautics.  Starting  in 
1979,  portions  of  the  famous  "Red  Book"  (The  Aerothermodynamics  of  Aircraft  Gas 
Turbine  Engines,  AFAPL-TR-78-52,  July  1978)  were  incorporated  into  the  Propulsion  II  and 
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Propulsion  Design  courses.  The  textbooks  used  in  all  three  courses  have  evolved  over  the 
past  decade  (see  Table  3)  and  continue  today. 

_ Tables  Textbooks  from  1976-85 _ 

Propulsion  I 

"Oates  Notes" 

Elements  of  Propulsion,  edited  by  Major  Jack  0  Mattingly* 

The  Jet  Age.  Forty  Years  of  Jet  Aviation,  edited  by  Boyne  and  Lopez* 

The  Aerothermodvnam  ICS  of  Gas  Turbine  and  Rocket  Propulsion,  by  Gordon  C,  Oates* 

Propulsion  II 

"Oates  Notes" 

Elements  of  Propulsion,  edited  bv  Maior  Jack  D  Mattingly* 

The  Aerothermodynamicsof  GasTurbineand  Rocket  Propulsion,  by  Gordon  C.  Oates* 

Aircraft  Gas  Turbine  Engine  Technology,  bv  Irwin  E  Traeoer* 

Red  Book  ("The  Aerothermodvnamics  of  Aircraft  Gas  Turbine  Engines,"  AFAPL-TR-78-52,  July  1978)** 

Propulsion  Design 

"Oates  Notes" 

Elements  of  Propulsion,  edited  by  Major  Jack  D.  Mattingly* 

The  Aerothermodynamicsof  Gas  Turbine  and  Rocket  Propulsion,  by  Gordon  C.  Oates* 

Aircraft  Gas  T urbine  Engine  Technology,  by  Irwin  E .  T raeger* 

Red  Book  ("The  Aerothermodynamicsof  Aircraft  Gas  Turbine  Engines,"  AFAPL-TR-78-52,  July  1978)** 
Fundamentals  of  Aircraft  Design  bv  Leland  M.  Nicolai** 

*  Current  text 

**  Reference 


The  cycle  analyses  of  "Oates  Notes"  were  incorporated  with  fundamental  gas 
dynamics  by  Major  Jack  Mattingly  in  1982  to  produce  Elements  of  Propulsion,  the  current 
text  of  the  Propulsion  I  course.  This  Academy  text  is  supplemented  with  readings 
about  the  early  development  of  gas  turbine  engines  from  The  Jet  Age.  Forty  Years  of 
Jet  Aviation,  edited  by  Boyne  and  Lopez  and  published  by  the  National  Air  and  Space 
Museum  Starting  this  past  spring,  the  Propulsion  I  course  began  use  of  The 
Aerothermodvnamics  of  Gas  Turbine  and  Rocket  Propulsion,  by  Gordon  C.  Oates  to 
supplement  the  material  contained  in  the  other  books. 

The  Propulsion  II  course  is  currently  using  both  Oates’  The  Aerothermodvnamics  of 
Gas  Turbine  and  Rocket  Propulsion  and  Aircraft  Gas  Turbine  Engine  Technology,  by  Irwin 
E.  Traeger.  The  text  by  Traeger  provides  excellent  material  on  component  and  engine 
construction  which  goes  hand-in-hand  with  our  cut-away  displays  of  J69  and  J85  engines. 
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The  Propulsion  Design  course  uses  numerous  handouts  that  are  supplemented  by 
Elements  of  Propulsion.  The  Aerothermodynamics  of  Aircraft  Gas  Turbine  Engines,  and 
Aircraft  Gas  Turbine  Engine  Technology.  The  methodology  of  the  constraint  analysis  and 
the  mission  analysis  used  in  the  Propulsion  Design  course  have  their  origins  in 
Fundamentals  of  Aircraft  Design  by  Leland  M.  Nicolai  (Ref.  3). 

IV.  The  Benefits  Cascade 

We  are  obviously  very  enthusiastic  about  the  U.S.  Air  Force  Academy  Department  of 
Aeronautics  propulsion  track,  and  especially  Propulsion  Design,  and  with  good  reason. 
Cadets  almost  invariably  rate  the  latter  "  the  best  course  they  have  ever  had."  The 
quantitative  results  have  always  looked  very  much  like  the  "real  world"  solution.  In  fact, 
the  students  often  find  useful  improvements  to  the  RFP  and/or  inevitable  characteristics 
of  the  resulting  engines  before  they  are  generally  known  elsewhere. 

These  activities  have  led  to  a  successful  summer  Propulsion  Workshop  at  the 
Academy,  which  is  a  blend  of  basic  cycle  training,  computer-aided  cycle  selection,  and 
informational  briefings  on  all  aspects  of  USAF  involvement  with  engine  research,  design, 
development,  acquisition,  operation,  and  maintenance.  Finally,  another  textbook  is  in 
the  works,  which  is  based  upon  and  will  capture  the  spirit  and  methodology  of 
Propulsion  Design.  This  textbook  is  programmed  for  publication  as  a  part  of  the  AlAA 
Education  Series  in  1986  as  Aircraft  Engine  Desiqn. 
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How  To  Fairly  Account  For  An  Excused  Exam 

Robert  C.  Winn* 

Abstract 

This  paper  describes  a  methodical  procedure  for  accurately  assigning  a  course 
grade  to  a  student  who  has  had  one  or  more  exams  excused.  Simply  excluding  that  exam 
from  the  student's  average  may  not  be  fair  to  that  student  or  to  the  other  students  in  the 
course,  depending  on  the  relative  difficulty  of  the  excused  exam.  An  equation  is 
presented  in  order  to  calculate  an  imaginary  score  for  the  excused  exam,  a  score  which  is 
consistent  with  the  student's  performance  on  past  exams.  Further,  a  computer  program 
called  FAKEGRD ,  which  automates  the  calculation  of  the  imaginary  scores  is  presented. 

I.  Introduction 

In  the  last  few  years,  the  workload  in  the  Department  of  Aeronautics  has  increased 
to  the  point  that  specific  efforts  to  reduce  it  have  been  considered.  One  policy  change 
which  was  implemented  was  the  elimination  of  the  requirement  for  make-up  exams. 
Historically,  make-up  exams  were  required  for  all  tests  and  course  quizzes  unless 
extraordinary  circumstances  existed.  This  put  a  large  burden  on  the  course  director 
because  the  additional  exams  were  time  consuming  to  prepare,  particularly  due  to  the 
difficulty  in  ensuring  equivalence  between  exams.  Eliminating  make-up  exams  clearly 
had  the  potential  to  reduce  instructor  workload. 

However,  before  excusing  a  cadet  from  a  graded  event,  the  instructor  must  answer 
a  very  important  question-How  does  he  or  she  assign  a  valid  course  grade?  Imbedded  in 
this  question  are  two  other  concerns.  First,  how  much  testing  is  enough  to  determine  a 
valid  grade  for  any  cadet,  and  second,  how  should  a  cadet's  grade  be  determined  if  an 
exam  is  excused?  For  the  purpose  of  this  paper,  I  will  assume  that  there  is  enough 
evaluation  of  each  student  to  allow  the  excusal  of  one  or  two  exams.  Determining  just 
how  much  testing  is  enough  to  establish  a  valid  course  grade  is  a  complex  question  and 
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will  not  be  addressed  here.  In  this  paper,  I  will  concentrate  on  the  question  of  how  to 
account  for  an  excused  exam  in  the  assignment  of  a  final  course  grade. 

II.  What  to  Do  About  an  Excused  Exam 


In  the  past,  when  a  student  was  excused  from  an  exam,  that  score  was  simply  left 
out  of  the  calculation  of  the  student  average.  If  all  exams  had  identical  means  and 
standard  deviations,  then  this  method  was  fine.  If,  however,  the  exams  varied  in 
difficulty,  then  some  real  problems  could  occur.  If  a  student  missed  a  very  difficult  exam, 
his  course  average  would  be  higher  than  it  should  be.  The  converse  was  also  true. 


Table  1 

Sample  Grades  with  Excused  Exams 


Student  A 

Student  B 

Course 

Average 

Course 

Std  Dev 

Maximum 

Possible 

Test#1 

70 

70 

70 

10 

100 

Test  #2 

— 

70 

60 

20 

100 

Test  #3 
AVG. 

70 

70.00 

70.00 

80 

70.00 

5 

100 

The  example  depicted  in  Table  1  shows  what  can  happen  if  missed  exams  are  not 
accounted  for  properly.  Both  students  have  exactly  the  same  average,  but  they  have  not 
performed  equally.  Test  #2  was  a  very  difficult  one  with  a  large  standard  deviation.  Test 
#3  was  very  easy  with  a  low  standard  deviation.  Student  A  did  average  on  Test  #1  but 
well  below  average  (2  standard  deviations  )  on  Test  #3.  Because  he  didn't  take  the 
difficult  Test  #2,  his  average  is  unrealistically  high.  Student  B  happened  to  miss  the  easy 
Test  #3.  As  a  result,  he  has  the  same  average  as  Student  A  even  though  he  equaled  or 
exceeded  the  mean  on  every  exam.  This  example  points  out  the  dangers  associated  with 
using  uncorrected  averages  when  an  exam  is  excused. 

If  you  examine  the  scores  in  the  above  example,  you  can  see  some  logic  that  will 
result  in  a  fairer  grade  for  both  students.  Student  A  was  a  total  of  10  points  below  the 
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I 

i 

I 

I  mean  on  the  two  exams  he  took.  More  accurately,  he  was  2  standard  deviations  below 

the  mean  for  the  two  exams  or  an  average  of  1  standard  deviation  below  the  mean  per 
exam.  Therefore,  the  best  estimate  of  what  his  performance  would  have  been  on  the 
exam  he  missed  is  1  standard  deviation  below  the  mean  on  the  exam  or  a  score  of  40%. 
Similarly,  Student  B  would  be  expected  to  earn  a  score  of  81.25%,  0.25  standard 
deviations  above  the  mean,  on  Test  #3.  By  using  these  imaginary  grades  in  the 
calculation  of  the  students'  averages.  Student  A  ends  up  with  an  average  of  60% 
compared  to  73.75%  for  Student  B.  The  effect  of  including  these  imaginary  scores  in  a 
student's  average  is  to  increase  the  weight  of  the  exams  he  did  take.  The  new  averages 
are  a  much  better  indication  of  student  performance  in  the  course. 

The  above  example  was  a  rather  drastic  one  because  of  the  large  differences  in 
standard  deviations  from  one  exam  to  another.  A  more  realistic  example  is  taken  from 
Fall  84  Aero  371  as  shown  in  Table  2. 


Table  2 

Grades  from  Fall  1984,  Aero  371 


Student  A 

Student  B 

Student  C 

Course 

Average 

Course 
Std  Dev 

Maximum 

Possible 

Test#1 

- 

48 

46 

42.6 

7  6 

50 

Test  #  2 

61 

59 

— 

64.4 

18.9 

100 

Test  #3 

86 

- 

70 

57.7 

15.4 

100 

Test  #4 

43 

49 

46 

43.5 

4.7 

50 

Test  #5 

66 

60 

52 

72.1 

14.5 

100 

AVG. 

73.  14% 

72.00% 

71.33% 

69.2% 

12.2 

The  same  basic  logic  used  in  the  earlier  example  can  be  used  in  this  more  complex  one. 
Student  A  receives  a  45  for  Test  #1  and  now  has  an  average  of  75.25%.  Student  B 
receives  a  57  for  Test  #3  and  now  has  an  average  of  68  25%.  Student  C  receives  a  64  for 
Test  #2  and  now  has  an  average  of  69.50%  These  imaginary  grades  corrected  the 
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misplacement  of  these  students  in  the  order  of  merit  which,  in  this  case,  has  76  students. 
Student  A  moved  from  29  to  26,  Student  B  moved  from  34  to  43,  and  Student  C  moved 
from  36  to  41.  Student  A  moved  up  because  he  missed  an  easy  exam.  Students  B  and  C 
moved  down  because  they  missed  hard  exams,  but  Student  C  moved  ahead  of  Student  B. 

Calculating  the  imaginary  grades  in  this  second  example  is  more  challenging  than 
in  the  first  example  because  not  all  exams  were  of  equal  weight.  To  handle  problems 
like  this  and  to  enable  automation  of  the  procedure,  I  developed  general  methodology 
for  the  calculation  of  the  imaginary  grades. 

III.  Methodical  Calculation  of  Imaginary  Grades 

You  can  calculate  the  imaginary  grades  quite  methodically  by  using  the  following 
equation: 


imag  mtsa  N 


n  =  l 


max,n 


N 


1 


X  ■  (x  -X  ) 
max,n  n  n 


a 

n 


where 

Ximag  =  the  imaginary  score  given  on  the  missed  exam 

^miss  =  mean  on  the  exam  which  was  missed 

Omiss  =  the  standard  deviation  on  the  exam  which  was  missed 

N  =  the  total  number  of  exams  taken  by  the  student  in  question 

^max.n  =  the  maximum  number  of  points  possible  on  Exam  #n 

Xn  =  the  score  achieved  by  the  student  in  question  on  Exam  #n 

Xn  =  the  mean  on  Exam  #n 

On  =  the  standard  deviation  on  Exam  #n 

If  you  carefully  study  this  equation,  you  can  get  some  meaning  from  eacn  term. 
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The  term  (xn  -  x^)  is  simply  the  difference  between  the  student's  score  and  the  mean  on 
Exam  #n.  By  dividing  this  term  by  the  standard  deviation,  you  get  the  number  of 
standard  deviations  away  from  the  mean  the  student  was  on  that  test.  Multiplying  this 
number  by  the  maximum  score  gives  a  number  which  is  a  weighted  difference  from  the 
mean.  By  adding  these  weighted  differences  for  each  test  that  was  taken,  you  arrive  at 
the  total  difference  from  the  mean  for  the  whole  course.  That  divided  by  the  maximum 
possible  points  for  the  exams  taken  yields  an  average  number  of  standard  deviations 
from  the  mean  for  this  student.  Multiply  that  number  by  the  standard  deviation  for  the 
missed  exam  and  add  the  result  to  the  mean  of  the  missed  exam,  and  you  arrive  at  the 
imaginary  score  for  that  exam. 

A  simplified  form  of  the  above  equation  results  if  all  exams  have  the  same 
maximum  number  of  points  available.  The  above  equation  reduces  to 

X  =X  .  X 

imag  miss  "  n  n  n 


This  equation  is  easier  to  understand  than  the  more  general  one.  Note  that  it  can  be 
used  on  the  first  example  in  this  paper,  but  not  on  the  second  example. 

I  wrote  a  computer  program  called  FAKEGRD  to  assist  in  the  implementation  of 
the  imaginary  grade  calculation.  A  listing  of  the  program  is  attached  as  Appendix 
1 .  This  program  reads  the  exam  data  from  the  course  data  disk,  performs  the  calculation 
of  the  imaginary  scores,  and  prints  the  results  for  each  student  who  missed  an  exam.  If 
the  imaginary  scores  look  reasonable,  they  should  then  be  included  on  the  course  data 


IV.  Cautions  on  the  Use  of  FAKEGRD 

I  have  a  few  recommendations  and  warnings  concerning  the  calculation  of 
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imaginary  grades  in  general. 

1.  The  imaginary  grades  should  be  determined  after  all  graded  work  except  the 
final  has  been  entered  to  the  course  data  disk.  At  the  course  director's  discretion, 
Instructor  Prerogative  (IP)  points  and  graded  reports  could  also  be  used  in  the  calculation 
of  the  imaginary  scores.  If  the  course  director  decides  not  to  use  all  of  the  grades  in  the 
course  file,  he  must  create  a  new  file  and  save  it  on  the  disk  under  a  new  name. 

2.  Do  not  give  an  imaginary  score  greater  than  100%  or  less  than  0%,  regardless 
what  the  FAKEGRD  program  tells  you.  By  calculating  an  imaginary  score  for  a  missed 
exam,  we  are  trying  to  estimate  what  the  student  would  have  scored  if  hehadtaken  the 
exam.  He  could  not  have  scored  greater  than  100%  if  he  had  taken  the  exam. 

3.  Imaginary  grades  can  and  should  be  calculated  at  mid-term,  but  these  grades 
should  only  be  used  qualitatively  to  help  determine  a  mid-term  letter  grade.  They 
should  not  be  stored  on  the  course  grade  disk  at  mid-term.  A  new  imaginary  grade  will 
be  calculated  at  the  end  of  the  semester  after  the  student  has  established  his  position  in 
the  class  standing. 

V.  Conclusion 

Faculty  workload  will  continue  to  be  a  problem.  As  we  take  steps  to  reduce  it  we 
must  be  careful  to  be  as  fair  as  possible  to  each  student.  The  imaginary  grade  calculation 
described  in  this  paper  will  ensure  that,  if  a  graded  event  is  excused,  a  fair  course  grade 
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PROGRAM  FAKEGRD,  28  AUG  85 


DIMENSION  XB(14),X(14,200),MISS(14),XM(14),SD(14),AM(14),STUD(200) 
CHARACTER  A*1 ,CRS*8,STUD»8,AM»1 ,EX»7,D«8,F»12,DAT»8 

SET  UP  THE  PROGRAM. 

WRITE  (*,1001) 

1001  FORMAT  (//////////,'  THIS  PROGRAM  IS  FOR  USE  IN  COURSES  WITH  LESS' 

1, '  THAN  200  STUDENTS.',/,'  IF  YOUR  COURSE  HAS  MORE  THAN  200', 

2'  STUDENTS,  USE  THE  "CORE"  VERSION  OF  FAKEGRD.',/, 

3'  IF  YOU  WANT  TO  EXIT  THIS', 

4'  PROGRAM  MOW,  TYPE  "STOP";  IF  NOT,  HIT  "RETURN".',/////////) 

READ  (*,2001)  A 

2001  FORMAT  (A1) 

IF  (A.EQ.'S')  STOP 
WRITE  (*,1002) 

1002  FORMAT  (/////,'  THIS  PROGRAM  WILL  PROVIDE  YOU  WITH  IMAGINARY', 

1'  SCORES  FOR  MISSED  EXAMS.',//,'  THE  RESULTS  WILL  BE  PRINTED,  SO', 
2'  TURN  ON  THE  PRINTER.',//,'  AFTER  YOU  GET  THE  RESULTS,  CHECK', 

3'  THEM  CAREFULLY.',/,'  IF  THEY  LOOK  GOOD,', 

4'  INPUT  THEM  MANUALLY  USING  THE  "GRADER"  PROGRAM.',/) 

WRITE  (*,1003) 

1003  FORMAT  ('  NOTE:  IN  DETERMINING  THE  IMAGINARY  GRADES,  YOU  SHOULD', 

r  USE  ALL',/,'  OF  THE  GRADES  (EXCEPT  THE  FINAL)  WHICH  HAVE' 

2, '  BEEN  QUANTITATIVELY',/,'  ESTABLISHED.  IP  POINTS  MAY  BE' 

3, '  USED  IF  YOU  WISH.  ',/,'  IF  YOU  DO  NOT  WANT  TO  USE  ALL', 

4'  OF  THE  GRADES,  CREATE  A  NEW',/,'  DATA  FILE  AND  RENAME', 

5'  IT.  IF  YOU  WANT  TO  LEAVE  THE',/,'  PROGRAM  NOW,  TYPE', 

6'  "STOP";  IF  NOT,  HIT  "RETURN". '//////) 

READ  (*,2001)  A 
IF  (A.EQ.'S')  STOP 
WRITE  (*,1004) 

1004  FORMAT  (///,'  TYPE  YOUR  COURSE  WITH  ".DAT"  APPENDED  (I.E.  ' 
r"AER0356.DAT").' ,/,'  NOTE:  IF  YOU  CREATED  A  NEW  DATA  FILE,' 

2'  TYPE  IN  THAT  NEW  FILE  NAME.',/////) 

READ  (*,2002)  F 

2002  FORMAT  (A  12) 

OPEN  (5,FILE=F) 

OPEN  (6,FILE='PRGRD' ,STATUS= ' NEW' ) 

:  READ  THE  PRELIMINARY  EXAM  DATA  FROM  FILE  "F". 

WRITE  (*,1005) 

1005  FORMAT  (/,'  NOW  PLEASE  INPUT  THE  NUMBER  OF  STUDENTS  IN  YOUR', 

1'  COURSE  AND  THE  NUMBER  OF  EXAMS  GIVEN.',////, 

2'  THE  NUMBERS  SHOWN  BELOW  ARE  THE  NUMBERS  YOU  SHOULD  TYPE  IN.',/) 
READ  (5,2003)  CRS.DAT 
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D  Line#  1  7  Microsoft  FORTRAN??  V3.10  05/03/83 

55  2003  FORMAT  (A8,/,8X,A8) 

56  93  WRITE  (‘.lOOe)  DAT 

5?  1006  FORMAT  (1X,A8,//'  TYPE  IN  THE  FIRST  NUMBER,  THE  NUMBER  OF', 

58  r  STUDENTS,  AND  HIT  RETURN.') 

59  READ  (»,2004)  NSTUDS 

60  2004  FORMAT  (14) 

61  WRITE  (»,100?) 

62  100?  FORMAT  ('  NOW  TYPE  IN  THE  SECOND  NUMBER,  THE  NUMBER  OF', 

63  1'  EXAMS,  AND  HIT  RETURN.') 

64  READ  (•,2004)  NEXAMS 

65  WRITE  (*,1008)  NSTUDS, NEXAMS 

66  1008  FORMAT  (///,'  DO  THESE  NUMBERS, ' ,14 , '  AND ',13,',  MATCH  THE', 

6?  1'  ONES  I  PRINTED?') 

68  READ  (*,2001)  A 

69  IF  (A.EQ.'Y')  GO  TO  95 

?0  GO  TO  93 

?1  95  IF  (NEXAMS.lt. 9)  GO  TO  100 

?2  WRITE  (*,1009) 

?3  1009  FORMAT  (//,'  YOU  SHOULD  HAVE  TYPED  "COMPRESS"  BEFORE  "FAKEGRD".' 

?4  1,/,'  IF  YOU  DID  NOT,  TYPE  "STOP",  THEN  "COMPRESS",  THEN  ', 

?5  2'"FAKEGRD".' ,/,'  IF  YOU  DID,  JUST  HIT  "RETURN".') 

?6  READ  (*,2001)  A 

??  IF  (A.EQ.'S')  STOP 

?8  C 

?9  C  READ  IN  ALL  OF  THE  EXAM  SCORES  FOR  EACH  STUDENT. 

80  C 

81  100  WRITE  (*,1010) 

82  1010  FORMAT  (//////////,'  THE  CALCULATIONS  ARE  UNDERWAY.',//, 

83  1'  SIT  BACK  AND  RELAX;  THIS  MAY  TAKE  A  WHILE.',/////////////) 

84  READ  (5,2005)  (XM(I) ,1=1 , NEXAMS) 

85  2005  FORMAT  (/,20(F3.0, IX)) 

86  TOT  =  0. 

8?  DO  150  1=1, NEXAMS 

1  88  150  TOT=TOT+XM(I) 

89  DO  200  N=1, NSTUDS 

1  90  200  READ  (5,2006)  STUD(N) , (X(I ,N) , 1=1 , NEXAMS) 

91  2006  FORMAT  (A8,21X,20F4.0) 

92  EX='  EXAM#' 

93  D=' - ' 

94  C 

95  C  CALCULATE  THE  MEANS  AND  STANDARD  DEVIATIONS 

96  C 


97 

DO  300  1=1 .NEXAMS 

1 

98 

STUDS  =  0. 

1 

99 

XB(I)  =  0. 

1 

100 

DO  250  N=1 .NSTUDS 

2 

101 

IF  (X(I,N).LT.O.)  GO  TO  250 

2 

102 

STUDS  =  STUDS  ♦  1. 

2 

103 

XB(I)  =  XB(I)  ♦  X(I,N) 

2 

104 

250 

CONTINUE 

1 

105 

300 

XB(I)  =  XB(I)/STUDS 

106 

DO  400  1=1 .NEXAMS 

1 

10? 

STUDS  =  0, 

1 

108 

SD(I)  =  0. 

77 
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109 

DO  350  N=1,NSTUDS 

2 

110 

IF  (X(I,N).LT.0.)  GO  TO  350 

2 

111 

STUDS  =  STUDS  +  1. 

2 

112 

SD(I)  =  SD(I)  +  (X(I,N)  -  XB(I))»»2 

2 

113 

350 

CONTINUE 

1 

114 

400 

SD(I)  =  SQRT(SD(I)/(STUDS  -  1.)) 

115 

C 

116 

C 

SEND  THE  HEADING  TO  THE  PRINT  FILE 

11? 

c 

1 18 

WRITE  (6,1011)  CRS 

119 

1011 

FORMAT  (////, 30X, »FAKE  GRADES  FOR  ',A8,//) 

120 

WRITE  (6,1012) 

121 

1012 

FORMAT  (•  NOTE:  ASTERISKS  (*)  IDENTIFY  THE  IMAGINARY', 

122 

1'  SCORES.',/) 

123 

WRITE  (6,1013)  (EX,I,I=1,NEXAMS) 

124 

1013 

FORMAT  (/,'  STUDENT  NEW  AVG' , 15(A?,I1 ) ) 

125 

WRITE  (6,1014)  (D,I=1,NEXAMS) 

126 

1014 

FORMAT  (' - ',15A8) 

12? 

C 

128 

C 

EVALUATE  THE  SUMMATIONS  NECESSARY  FOR  THE  CALCULATION. 

129 

C 

130 

DO  500  N=1,NSTUDS 

1 

131 

TERM=0 

1 

132 

XMAX=0 

1 

133 

N0MISS=0 

1 

134 

K  =  1 

1 

135 

DO  420  I=1,NEXAMS 

2 

136 

AM(I)='  ' 

2 

13? 

IF  (X(I,N).GE.O.)  GO  TO  410 

2 

138 

MISS(K)=I 

2 

139 

NOMISS=l 

2 

140 

K=K+1 

2 

141 

GO  TO  420 

2 

142 

410 

XMAX=XMAX+XM(I) 

2 

143 

TERM=TERM+XM(I)*(X(I,N)-XB(I))/SD(I) 

2 

144 

420 

CONTINUE 

1 

145 

IF  (NOMISS.EQ.O)  GO  TO  500 

1 

146 

K=K-1 

1 

14? 

C 

1 

148 

C 

CALCULATE  THE  IMAGINARY  GRADE(S). 

1 

149 

c 

1 

150 

DO  440  1=1, K 

2 

151 

M=MISS(I) 

2 

152 

X(M,N)=XB(M)+SD(M)»TERM/XMAX 

2 

153 

AM(M)='*' 

2 

154 

440 

CONTINUE 

1 

155 

c 

1 

156 

c 

SEND  THE  CALCULATED  INFORMATION  TO  THE  PRINT  FILE. 

1 

15? 

c 

1 

158 

AVG=0. 

1 

159 

DO  450  I=1,NEXAMS 

2 

160 

450 

AVG=AVG+X(I,N) 

1 

161 

AVG=100.»AVG/T0T 

1 

162 

WRITE  (6,1015)  STUD(N),AVG,(X(I,N),AM(I),I=1,NEXAMS) 

05:10:13 

Microsoft  FORTRAN??  V3.10  05/03/83 


D  Line#  1  ? 

1  163  1015  FORMAT  ( IX, A8,F? .2, 10(F?.0,A1 ) ) 

1  164  500  CONTINUE 

165  WRITE  (6,1014)  (D,I=1 .NEXAMS) 

166  WRITE  (6,1016)  (XM(I) ,1=1 , NEXAMS) 

16?  1016  FORMAT  ('  MAX’ , 1 3X, 1 4F8. 1 ) 

168  WRITE  (6,101?)  (XB(I) ,1=1 , NEXAMS) 

169  101?  FORMAT  ('  MEAN ' , 1 2X, 14F8.2) 

1?0  WRITE  (6,1018)  (SD(I), 1=1, NEXAMS) 

1?1  1018  FORMAT  ('  STD  DEV ,9X, 14F8.2) 

1?2  WRITE  (6,1019) 

1?3  1019  FORMAT  (//,'  ENSURE  ACCURACY  BY:',/, 

1?4  1'  1.  LOOKING  FOR  "WAY-OUT"  NUMBERS  LIKE  ZERO.',/, 

1?5  2'  2.  CHECKING  THE  FAKE  GRADES  FOR  REASONABLENESS.',/, 

1?6  3'  3.  ENSURING  THAT  THERE  IS  A  FAKE  GRADE  FOR  EVERY' 

1??  4'  MISSING  GRADE.',/, 

1?8  5'  4.  COMPARING  THE  MEAN  AND  STD  DEV  HERE  AGAINST  A', 

1?9  6'  "GRADER"  OUTPUT.',//) 

180  WRITE  (6,1020) 

181  1020  FORMAT(/,'  NOTE:  DO  NOT  INPUT  FAKE  GRADES', 

182  1'  TO  THE  "GRADER"  PROGRAM  AT  THE  PROG.',/,?X, 

183  2'  JUST  USE  THIS  INFORMATION  TO  HELP  YOU  ASSIGN  PROG  GRADES.',/ 

184  3,?X,'  INPUT  THE  FAKE  GRADES  USING  "GRADER"  AT  THE  END  OF  THE’, 

185  4'  SEMESTER  ONLY.') 

186  CLOSE  (5) 

18?  CLOSE  (6) 

188  WRITE  (*,1021) 

189  1021  FORMAT  (/,'  I  AM  SENDING  THE  RESULTS  TO  THE  PRINTER  NOW.', 

190  1///,'  IT  HAS  BEEN  A  PLEASURE  WORKING  WITH  YOU.',/////////) 

191  STOP 

192  END 


Name 

Type 

Offset  P  Class 

A 

CHAR»1 

13360 

AM 

CHAR»1 

13026 

AVG 

REAL 

15418 

CRS 

CHAR»8 

14486 

D 

CHAR»8 

15151 

DAT 

CHAR»8 

14494 

EX 

CHAR»? 

15144 

F 

CHAR* 12 

142?4 

I 

INTEGER*4 

15088 

K 

INTEGER*4 

15402 

M 

INTEGER»4 

15414 

MISS 

INTEGER»4 

11314 

N 

INTEGER»4 

15120 

NEXAMS 

INTEGER*4 

14?04 

NOMISS 

INTEGER*4 

15398 

NSTUDS 

INTEGER»4 

14606 

SD 

REAL 

113?0 

SORT 

INTRINSIC 

STUD 

CHAR«8 

11426 

STUDS 

REAL 

15164 

TERM 

REAL 

15390 

D  Line! 

1  ? 

TOT 

REAL 

15112 

X 

REAL 

114 

XB 

REAL 

2 

XM 

REAL 

58 

XMAX 

REAL 

15394 

Name 

Type 

Size 
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PROGRAM 


Pass  One  No  Errors  Detected 
192  Source  Lines 
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